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ABSTRACT 

This is the third paper in a series that reports on our investigation of the clustering properties of 
AGNs identified in the i? 054 T All- Sky Survey (RASS) and Sloan Digital Sky Survey (SDSS). In this 
paper, we extend the redshift range to 0.07 < z < 0.50 and measure the clustering amplitudes of both 
X-ray and optically-selected SDSS broad-line AGNs with and without radio detections as well as for 
X-ray selected narrow-line RASS/SDSS AGNs. We measure the clustering amplitude through cross- 
correlation functions (CCFs) with SDSS galaxies and derive the bias by applying a halo occupation 
distribution (HOD) model directly to the CCFs. We find no statistically convincing difference in the 
clustering of X-ray and optically-selected broad-line AGNs, as well as with samples in which radio- 
detected AGNs are excluded. This is in contrast to low redshift optically-selected narrow-line AGNs, 
where radio-loud AGNs are found in more massive halos than optical AGNs without a radio-detection. 
The typical dark matter halo masses of our broad-line AGNs are log {MY^M\il[h~^ -^o]) ~ 12-4 — 13.4, 
consistent with the halo mass range of typical non-AGN galaxies at low redshifts. We find no significant 
difference between the clustering of X-ray selected narrow-line AGNs and broad-line AGNs. We 
confirm the weak dependence of the clustering strength on AGN X-ray luminosity at a ~2(T level. 
Finally, we summarize the current picture of AGN clustering to z ~ 1.5 based on three dimensional 
clustering measurements. 

Subject headings: galaxies: active ~ X-rays: galaxies - cosmology: large-scale structure of Universe 



1. INTRODUCTION 

Galaxy clustering measurements of large area sur- 
veys, such as the 2dF G alaxy Redshift Survey (2dF- 
GRS, I Colless et al.) [200lh. S loan Digital Sky Survey 
(SDSS. Abaz aiian et al.ll2009[ ). Deep Extrag alactic Evo- 
lutionary Probe 2 Redshift Survey (DEEP2. IDavis et al.l 
2003), AGN a nd Ga laxy Evolution Survey (AGES, 
Kochanek et aP l200l. and VIMOS-VLT Deep Sur- 
vev (VVDS. iLe Fevre et al.l [2005.1 with many thou- 
sands of galaxies have quantified the clustering de- 
pendence on galaxy properties such as morphologi- 
cal type, l uminosity, spectral type, and redshift to 
z ^ 1 (e.g., jtt jorberg et al."2002': 'MadEwick et al.' '2003'; 
Zehavi et alT[2005b : Mcncux ct al. 2006; Coil ct al. 2008; 
Meneux et al.ll2009l) . The general trend is that red, more 
massive, brighter, and/or elliptical galaxies cluster more 
strongly than blue, less massive, fainter, and/or spiral 
galaxies. Halo occupation distribution (HOD) modelling 
of these results has shown that there is an almost linear 
increase in the mean number of satellite (non-central) 
galaxies as a fun ction of increasing dark matter hal o 
(DMH) mass (e.g.. lZheng et al.|[2007t[Zehavi et al.|[20TTI ). 

Compared to projected angular clustering measure- 
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ments, three dimensional (3D) clustering measurements 
have greater statistical accuracy and minimize system- 
atic errors associated with m odel assumptions in using 
Limber's de-projection ([Limb er 1954). However, 3D cor- 
relation functions require extremely large numbers of ob- 
jects with secure redshift information. 

Clustering measurements of the auto-correlation func- 
tion (AGE) of active galactic nuclei (AGNs) reveal 
how AGNs are spatially distributed in the Universe. 
These measurements can constrain the AGN host 
galaxies, determine the typical DMH mass in which 
AGNs reside and the distribution of AGNs with DHM 
mass, test theoretical model predictions, and address 
which ph ysical processe s are triggering AGN activ- 
itv (e.g.. (Porciani et all [200l iGilli et al.l [20051 
l Yang et all I200fif ICoil eTalT |2007l: ICoil et al.l 



[Ross ct al. 2009'; Krumoc ct al. 2010a; Cappe lluti et al.l 
2010; Mivaj i et al., 2011; AUevato et aL 2011). However, 
a major challenge in AGN clustering measurements is to 
overcome the limitation caused by the low space density 
of AGNs. 

Narrow emission lin e di agnostics (e.g., 
iVeilleux fc Osterbrockl Il987t [Kauffmann et al. 200|) 
have been used to identify narrow-line AGNs at low 
redshifts. Large sky area surveys with extensive spec- 
troscopic follow-up programs have recently allowed 
narrow-line AGN clust ering measurem ents for samples 
of up to 90 , 000 A G Ns (iLi et al.l[2006[). Studies such as 
'Wak e et al.l (|20p). 'Li et al.' (2006*), 'Magliocchetti et al.| 
(2004f^ and HMandelbaum et al, (,2009l) explore the 
clustering strength of low redshift {z < 0.3) narrow-line 
AGNs samples with respect to various AGN properties 
(e.g., black hole mass, radio emission). Optically- 
selected narrow-line AGNs have a clustering strength 



2 



Krumpe et al. 



similar to all galaxies (both blue and red). Narrow- line 
radio-loud AGNs cluster more strongly than AGN 
without radio-emission and identical samples of qui- 
escent galaxies. Consequently, denser environment of 
galaxies increase the probability of hosting a narrow-line 
radio-loud AGN. 

Broad-line AGNs are, in general, more luminous than 
narrow-line AGNs. Moreover, their energy production 
dominates over the host galaxy star light. Large area 
optical surveys have spectroscopically identifie d tens 
of th ousands of broad-line AGNs (e.g., Schneider et al.l 
I2OII rCroom et aLH200l . The co-moving number den- 
sity of broad-line AGNs is very low in the low redshift 
universe {z < 0.5). The situation improves at higher red- 
shifts as the co-moving AGN number distribution peaks 
at 2 ^ 2. Furthermore, a given observed area corre- 
sponds at higher redshift to a larger observed co-moving 
volume. These two facts favor bro ad-line AGN cluster- 
ing measurements at z > 0.5 (e.g.. jPorciani et al.ll2004l : 
iCroom et all [20051: iRoss et aLl^OOl . 

X-ray surveys have also been used to select AGNs. 
These surveys cover much less sky area, e.g., ~0.1-10 
deg^, than optical surveys. Further, they sample lower 
AGN luminosities than optical surveys and contain a sig- 
nificant fraction of obscured AGNs. In order to probe 
the required large observed co-moving volume needed 
for clustering measurements using these relatively small 
sky areas, 3D clustering measurements of X-ray selected 
AGNs have mainly been conducted at z > 0.8 (e.g . , 
IGilh et al.ll2005L[2009l:lYang et al.l[200l ICoil et al.l[2007h . 

Optically-selected AGNs at redshifts of z ^ 1.0 — 
1.5 appear to reside in somewhat low er host DMH 
mass e s (Mdmh 'y 10^^ " -^^ Mf^, e.g., JPorciani et all 
[200I ICoil et al.1 [20071: IRoss et al.l [2009h than XMM- 
Newton and Chandra selected AGN samples (Mdmh ~ 
-^013-13.5^-1 lYang et al.l [2OO6I: IGilh ef^ 

[200a [2009t [Coil et al.ll2009HAllevato et al.ll2011l) . Possi- 
ble explanations of the differences in the clustering sig- 
nals include either the presence of a large fraction of X- 
ray absorbed narrow-line type II AGNs in the XMM- 
Newton and Chandra samples or the different luminosi- 
ties of the AGN samples. At lower redshifts 3D broad- 
line X-ray AGN clustering measurements have been asso- 
ciate d with larg e uncertainties, due to small sample sizes 
(e.g.. lMullis et aL 2004). 

To achieve small er uncerta in ties i n AGN cluster- 
ing measurements, ICoil et al.l ([20091) use the cross- 
correlation function (CCF) of AGNs with a tracer set 
of galaxies that contains a large number of objects. The 
AGN ACF is then inferred from the CCF, which has 
many more pairs at a given separation and hence signifi- 
cantly reduces the uncertainties in the spatial correlation 
function compared with the direct measurements of the 
AG N ACF. 

In iKrumpe et al.l (|2010a| ) (hereafter paper I) we use 
the cross-correlation technique to calculate the CCF be- 
tween broad-line X-ray selected AGNs from the ROSAT 
AU Sky Survey (RASS) and SDSS Luminous Red Galax- 
ies (LRGs) in the redshift range 0.16 < z < 0.36. The 
potential of RASS, which is currently the most sensitive 
X-ray all sky survey, can only be maximally exploited 
when it is combined with other large-area surveys such 
as SDSS. The unprecedented low uncertainties of the in- 



ferred broad-line AGN ACF from the RASS/SDSS com- 
bination allows us to split our sample into subsamples ac- 
cording to their X-ray luminosities. From this work, for 
the first time, we report an X-ray luminosity dependence 
of broad-line AGN clustering in that higher luminosity 
AGNs cluster more strongly than their lower luminosity 
counterparts. We conclude that low luminosity broad- 
line RASS/SDSS AGNs cluster similarly to blue galaxies 
at the same redshift, while high luminosity RASS/SDSS 
AGNs cluster similarly to red galaxies^ 

In our second paper (Miyaii et al. 20111 : hereafter pa- 
per II), we apply a halo occupation distribution (HOD) 
modeling technique to the AGN-LRG CCF in order to 
move beyond determining the typical DMH mass based 
on the clustering signal strength and instead constrain 
the full distribution of AGNs as a function of DMH mass. 
To do this, we develop a novel method of applying the 
HOD model directly to the CCF. The HOD modeling 
significantly improves the analysis of the CCF because it 
properly uses the Fourier transformed linear power spec- 
trum in the "two halo term" as well as the non-linear 
growth of matter in the "one halo term" through the for- 
mation and growth of DMHs. This results in significant 
improvements over the standard method, which is used in 
paper I, of fitting both regimes with a phenomenological 
power law. One of the important results of this analysis 
is that at 0.16 < z < 0.36 the mean number of satellite 
AGNs in a DMH does not proportionally increase with 
halo mass, as is found for satellite galaxies. The AGN 
fraction among satellite galaxies actually decreases with 
increasing DMH mass beyond Mdmh 10^^ Mq. 

In this paper we extend the scope of our previous pa- 
pers to both somewhat lower and higher redshifts to 
obtain broad- line AGN clustering results at z < 0.5, 
where very precise narrow-line AGN clustering measure- 
ments exist but broad-line AGN clustering is poorly con- 
strained. This is crucial for studying the possible evolu- 
tion and luminosity dependence of broad-line AGN clus- 
tering from low to high redshifts. The dominant process 
that triggers AGN activity could be a function of red- 
shift and/or halo mass, which may be reflected in the 
clustering properties. We also study the clustering sig- 
nal of both X-ray and optically-selected broad-line AGN 
samples, and test whether the exclusion of radio-detected 
broad-line AGNs changes our results. As the same statis- 
tical method and galaxy tracer sets are used to infer the 
clustering signal for X-ray and optically-selected broad- 
line AGNs, we can explore differences in the clustering 
among the different selection techniques with low system- 
atic uncertainties. Furthermore, we derive bias parame- 
ters by applying the HOD modeling directly to all CCFs. 
In Paper II we show that using power law fit results even 
in the non-linear regime, as is commonly done in the liter- 
ature, is appropriate to detect difference in the clustering 
properties between different samples. However, the de- 
rived bias parameters and DMH masses based on these 
fits should be interpreted with caution as the fit does not 
only consider the linear regime (two halo term), but also 
the non-linear regime (one halo term). Consequently, 
here we derive the bias for each AGN sample using HOD 
modeling of the CCF. Full detailed results of the HOD 
modeling of the CCFs presented in this paper will be 
given in a future paper (Miyaji et al. in preparation). 

This paper is organized as follows. In Section 2 we 
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describe the properties of the different galaxy tracers 
sets used at different redshifts, while Section 3 gives 
the details of our different AGN samples. In Section 
4 we briefly summarize the cross-correlation technique, 
how the AGN ACF is inferred from this, and present 
our results. We apply the HOD modeling in Section 5. 
Our results are discussed in Section 6 and we conclude 
in Section 7. Throughout the paper, all distances are 
measured in co-moving coordinates and given in units 
of Mpc, where h = i?o/100kms~^ Mpc~^, unless 
otherwise stated. We use a cosmology of rim — 0.3, 
r^A = 0.7, and (Ts{z = 0) = 0.8, which is consis tent with 
the WMAP data release 7 ([Larson et al.ll2011l: Table 3). 
The same cosmology is used in papers I & II. Luminosi- 
ties and absolute magnitudes are calculated for h = 0.7. 
We use AB magnitudes throughout the paper. All un- 
certainties represent to a la (68.3%) confidence interval 
unless otherwise stated. 

2. GALAXY TRACER SETS 

The crucial ingredient in the cross-correlation method 
is the tracer set, a sample with a large number of ob- 
jects that traces the underlying dark matter density dis- 
tribution. The properties of the tracer set determines 
the redshift range over which the method can be ap- 
plied. The AGN samples of interest are necessarily lim- 
ited to the same redshift range and geometry as the corre- 
sponding tracer set. A s the RASS/SDSS selected AGNs 
([Anderson et al.l 120071 ) are based on the SDSS data re- 
lease 5 (DR5), we consequently limit the tracer sets to 
the same survey geometry when we compute the clus- 
tering measurements of the X-ray selected AGN. The 
SDSS geometry and c ompleteness are expressed in terms 
of spherical polygons fHam ilton fc Tegmarkll200l . This 
file is not publicly-available for DR5, therefore, we use 
the latest version available prior to DR5: the DR4-I- ge- 
ometry filj^. Howeve r, optically-selected SDSS AGNs 
([Schneider et al.ll2010D make use of the full SDSS survey 
(DR7). Consequently, we consider tracer sets from the 
DR7 geometr}|3 whenever we compute CCFs of optically- 
selected SDSS AGNs. 

In the redshift range of 0.07 < z < 0.16 we use 
SDSS main galaxies for the tracer set, while SPSS lu- 
minous red galaxies (LRGs, lEisenstein et al.l [20011 1 are 
used for the corresponding cross-correlation sample at 
0.16 < z < 0.36 (same as for paper I). Very luminous red 
galaxies are used as a tracer set at (3.36 < z < 0.50. We 
will refer to the latter sample as the 'extended LRG sam- 
ple'. Above z ~ 0.5 the number of galaxies with spec- 
troscopic redshifts in SDSS decreases dramatically and 
does not allow the selection of further tracer sets with a 
sufficient density of objects. In the following subsections 
we describe in detail the extraction of the various tracer 
sets and how we account for SDSS fiber collisions. 

2.1. SDSS Main Galaxy Sample 

The SDSS Main Galaxy sample ([Strauss et al.ll2002t ) 
is drawn fro m the NYU Value-Added Galaxy catalog 
(NYU VAGC. lBlanton et al.l l2( 305aHPadmanabha,n et at 
mM ). based on the SDSS DR7 ([Abazaiian et al.ll2009i r 

® http : //sdss ■ physics ■ nyu ■ edu/ lss/dr4plus | 

http: //sdss .physics .nyu. edu/lss/dr72J 
* http: //sdss .physics .nyu. edu/vagc/ 



The photometric data covers an area of 10417 deg^, 
while the spectroscopic data covers 7966 deg^. Absolute 
magnitudes, b ased on P etrosian fiux es, are j'T -corrected 
to z = 0.1 (Bla nton et al . 2003; Blant on fc Row eis 2007), 
which is close to the median red shift o f our sample. We 
follow a scheme similar to Z ehavi et al.l (j2005b. ') , who use 
the NYU VAGC to measure the clustering of various 
luminosity and color-selected galaxy subsamples. Sim- 
ilarly, we limit our sample to 14.5 < r < 17.5. The 
bright limit avoids incompleteness due to galaxy deblend- 
ing, and the faint limit accounts for the slightly- varying 
r-band magnitude limit over the SDSS area (nominal 
value r ^ 17.7). The restriction of r < 17.5 ensures 
a uniform flux limit throughout the whole SDSS sur- 
vey. In addition, we create a volume-limited galaxy sam- 
ple by selecting objects with an absolute magnitude of 
-22.1 < i < -21.1. Finally we limit the redshift 
range to 0.07 < z < 0.16. 

Applying the above-mentioned selection criteria and 
considering only SDSS DR7 areas with a spectroscopic 
completeness ratio of >0.8, we select 68273 galaxies 
with spectroscopic redshifts from the corresponding NYU 
VAGC. The DR7 area covered by restricting the spec- 
troscopic completeness ratio to >0.8 is 7670 deg^. The 
properties of this DR7 SDSS main galaxy sample are 
summarized in Table |T] 

As described above, the X-ray selected RASS/SDSS 
AGN samples are based on DR5. Therefore, we further 
reconflgure the DR7 tracer sets (in this case the SDSS 
main galaxy sample) to the DR4-I- geometry to deflne 
a common survey geometry to use when measuring the 
clustering of X-ray selected AGNs. The restriction to 
the DR4-I- survey area with a DR7 spectroscopic com- 
pleteness ratio of >0.8 corresponds to an area of 5468 
deg2. Table H] lists the properties of the DR4-h SDSS 
Main Galaxy sample. 

2.1.1. Accounting for the SDSS Fiber Collision 

An operational constraint of the SDSS spectroscopic 
program is that two fibers cannot be placed closer than 
55 arcsec on a single plate. Overlapping spectroscopic 
plates compensate partially for the effect. However, ~7% 
of the target galaxies cannot be spectroscopically ob- 
served because of fiber collisions. This observational bias 
is corrected by assigning to each galaxy that has not been 
observed the redshift of the ir nearest neighbor w ith a 
spectroscopic SDSS redshift ([Blanton et al.ll2"005aD . 

Although one might be concerned that this simple 
method could overcorrect and result in too many close 
galaxy pairs at the same redshi ft, which would then 
distort clustering measurements, iZehavi et al.l ([2005bl ) 
demonstrate that this correction procedure works very 
well. They use ACDM A^-body simulations and design 
three galaxy samples and measure the correlation func- 
tion for three samples: i) from the full simulated galaxy 
distribution, ii) from simulated SDSS data including fiber 
collision losses and not correcting for it, and iii) from 
simulated SDSS data that corrects for the fiber collision 
by assigning the redshift of their nearest spectroscopic 
neighbor. They verify that the differences between i) 
and iii) are much smaller than the statistical uncertain- 
ties down to scales of Vp ~ 0.1 Mpc, while ii) under- 
estimates the correlation function at scales rp < 1 
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Mpc. Therefore, we use the same fiber correction pro- 
cedure for our main galaxy sample. The NYU VAGC0 
provides this information through an SDSS fiber collision 
corrected galaxy sample. 

2.1.2. Construction of the Random Main Galaxy Sample 

The random sample is another crucial ingredient re- 
quired for measuring the correlation function. The pur- 
pose is to create a randomly-distributed sample of ob- 
jects that exactly matches all observational biases (win- 
dow function, redshift distribution, etc.) of the observed 
sample. We follow the procedure of paper I (for details 
see paper I, Section 3.1) and generate a set of random 
RA and Dec values within DR7 areas with a spectro- 
scopic completeness ratio of >0.8, populate areas with 
higher completeness ratios more than areas with lower 
completeness ratios, smooth the obs erved redshift dis- 
tribu tion by applying a least-square (jSavitzkv fc Golavl 
[1961 low -pass filter, and use this smoothed redshift pro- 
file to randomly assign redshifts to the objects in the 
sample. 

The number of objects in a random catalog is chosen 
to have an adequate number of pairs in the CCF at the 
smallest scales measured here. For clustering measure- 
ments with the main galaxy sample the random cata- 
log contains 100 times as many objects as the observed 
sample. The random catalog of the DR4+ SDSS main 
galaxy sample also contains 100 times more objects than 
the data and follows the same exact procedure described 
for the DR7 geometry, except that we restrict that area 
to the DR4-I- geometry. 

2.2. SDSS Luminous Red Galaxy Sample 

The selection of the SDSS LRG sample is described 
in detail in Section 2.1 of paper I. Here we briefly sum- 
marize the sample selection. We extract LRGs from the 
web-based SDSS Catalog Archive Server Jobs SvstenFl 
using the flag 'galaxy jed', which is based o n the selec- 
tion criteria defined in lEisenstein et aTl ()2001f ) . We verify 
that the extracted objects meet all LRG selection criteria 
and create a volume-limited sample with 0.16 < z < 0.36 
and -23.2 < M°-3 < -21.2, where M^-^ is based on 
the extinction-corrected r*^^.^^ magnitude to construct 
the /c-corrected and passively evolved rest-frame ffpetro 
magnitudes at z = 0.3. We consider only LRGs that fall 
into the SDSS area with a DR7 spectroscopic complete- 
ness ratio of >0.8 and have a redshift confidence level of 
>0.95. 

The correction for the SDSS fiber collision in the SDSS 
LRG sample is slightly different from that for the SDSS 
main galaxy sample. LRGs exhibit very well-defined 
spectra dominated by an old stellar population that 
evolves very slowly. The reduced scatter in the spectral 
energy distribution (SED) of LRGs results in much lower 
photometric redshift uncertainties than the estimates for 
main galaxies, which can have ongoing star formation 
and therefore have a wider distribution of SEDs. We 
make use of the precise LRG photometric redshifts to 
correct for the SDSS fiber collision in the following man- 
ner. We select from the SDSS archive all LRGs that pass 
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Fig. 1. — Absolute magnitude versus redshift for the extended 
SDSS luminous red galaxy sample (0.36 < z < 0.50). The absolute 
magnitude is based on the extinction-corrected r*^^.^^ magnitude, 
passively evolved to 2 = 0.3. The different selection criteria for low 
and high redsh ift LRGs are visible in the upper right corner (cut I 
and cut II, see lEisenstein et al.ll200lD . 

the pure photometric-based LRG selection criteria. We 
identify photometric LRGs that are closer than 55 arcsec 
to a spectroscopic observed LRG in our redshift and ab- 
solute magnitude range. We then assign a redshift using 
the following steps: we accept the spectroscopic redshift 
of the LRG even if its redshift confidence level is <0.95. If 
there is no spectroscopic redshift available for the object, 
we give the photometric LRG the same redshift as the 
spectroscopic neighbor LRG (within a 55 arcsec radius) 
if 

I ■2'spcc,j ^photo,i 1^ f^^photo,i,lcr- (1) 

If Equation [T] is not fulfilled, we assume that the pho- 
tometric redshift of the LRG is the correct redshift. 
A redshift is assigned only if the object meets the se- 
lection criteria to construct a volume-limited sample: 
0.16 < z < 0.36 and -23.2 < M°-^ < -21.2. Approx- 
imately 2% of the all LRGs in our sample are assigned 
redshifts. The properties of the sample are shown in Ta- 
ble [2 

The construction of the random catalogs is identical 
to the procedure described in Section 12.1 21 (for details 
see paper I, Section 3.1), except that the LRG random 
catalogs contain 200 times as many objects as the real 
DR7 and DR4+ LRG samples. This is a compromise 
between the required computation time to calculate the 
correlation functions and having sufficient counts on the 
smallest scales to avoid introducing noise. More objects 
are required at higher redshift for a given sky area to ac- 
count for the fact that with increasing redshift the same 
angular distance on the sky corresponds to larger physi- 
cal co-moving separations. 

2.3. Extended SDSS Luminous Red Galaxy Sample 

In order to extend our clustering measurements to 
higher redshifts, we create an "extended SDSS LRG sam- 
ple" over the redshift range of 0.36 < z < 0.50. The 
extraction of the sample and the correction for the SDSS 
fiber collision problem is identical to the SDSS LRG sam- 
ple (see Section [ 



ct.nMeiiiiyp4tiis.z^feiMpJet^TOuld also be volume- limited. 
However, that would require an absolute magnitude cut 
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TABLE 1 

Properties of the SDSS Galaxy Tracer Sets and the AGN Samples. 



Sample 


SDSS 




M, logLx range 


Sample 


(n) 


(M, logLx) 


Name 


Geometry 


2:-range 


(mag, ergs~^) 


Size 


{h^ Mpc-3) 


(2) (mag, erg s"-') 



SDSS Tracer Sets 



Main galaxy 


DR7 


0.07 < 


z 


< 0.16 


-22.1 < A/01 ^ _2i 1 


68273 


9.8 X 10- 


.4 


0.13 


-21.41 


sample 


DR4-f 


0.07 < 


z 


< 0.16 


-22.1 < AfO i < -21.1 


48994 


9.8 X 10" 


.4 


0.13 


-21.41 


LRG sample 


DR7 


0.16 < 


z 


< 0.36 


-23.2 < AfO-^ < -21.2 


65802 


9.8 X 10" 


. 5 


0.28 


-21.71 




DR4-I- 


0.16 < 


z 


< 0.36 


-23.2 < A/0-3 < -21.2 


45899 


9.6 X 10- 


- 5 


0.28 


-21.71 


Extended LRG 


DR7 


0.36 < 


z 


< 0.50 


-23.2 < A/0-3 ^ _2i.7 


28781 


3.9 X 10- 


- 5 


0.42 


-22.04 


sample 


DR4-I- 


0.36 < 


z 


< 0.50 


-23.2 < A/J'-^ < -21.7 


19831 


3.8 X 10- 


-5 


0.42 


-22.04 






X-ray Selected AGNs - RASS/SDSS AGN 














total RASS-AGN 


DR4-I- 


0.07 < 


z 


< 0.16 


43.05 < logLx < 44.12 


629 


5.2 X 


10- 


-5 


0.12 


43.59 


total RASS-AGN (rq) 


DR4-I- 


0.07 < 


z 


< 0.16 


43.04 < logLx < 44.04 


504 


4.5 X 


10- 


-5 


0.12 


43.55 


low Lx RASS-AGN 


DR4-I- 


0.07 < 


z 


< 0.16 


42.92 < logLx < 43.54 


293 


4.5 X 


10- 


-5 


0.11 


43.25 


high Lx RASS-AGN 


DR4-I- 


0.07 < 


z 


< 0.16 


43.54 < logLx < 44.27 


336 


7.0 X 


10- 


-6 


0.13 


43.89 


low Lx RASS-AGN (rq) 


DR4-I- 


0.07 < 


z 


< 0.16 


42.93 < logLx < 43.54 


253 


4.0 X 


10- 


-5 


0.11 


43.26 


high Lx RASS-AGN (rq) 


DR4-I- 


0.07 < 


z 


< 0.16 


43.54 < logLx < 44.16 


251 


5.2 X 


10- 


-6 


0.13 


43.85 


narrow linp R ASS- AON 


DR4-I- 


0.07 < 


z 


< 0.16 


42.81 < logLx < 43.81 


194 


6.5 X 


10- 


-4 


0.11 


43.32 






0.16 < 
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< 0.36 


43.69 < logLx < 44.68 




n V 
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1 n- 

lU 
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n 9^ 
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z 
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loo ( 


4.9 X 


10- 


-5 


U.ZO 


44. io 


low Lx RASS-AGN 


DR4-I- 


0.16 < 


z 


< 0.36 


43.62 < logLx < 44.29 


990 


5.8 X 


10- 


-5 


0.24 


43.94 


high Lx RASS-AGN 


DR4-I- 


0.16 < 


z 


< 0.36 


44.29 < logLx < 44.87 


562 


1.2 X 


10- 


-6 


0.28 


44.58 


low Lx RASS-AGN (rq) 


DR4-I- 


0.16 < 


z 


< 0.36 


43.64 < logLx < 44.29 


883 


4.8 X 


10- 


-5 


0.24 


43.95 


high Lx RASS-AGN (rq) 


DR4-I- 


0.16 < 


z 


< 0.36 


44.29 < logLx < 44.82 


454 


1.0 X 


10- 


-6 


0.28 


44.55 


nowi-Mir linn T? A A 1\T 
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z 


< 0.36 
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■6 
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total RASS-AGN 


DR4-I- 


0.36 < 


z 


< 0.50 


44.25 < logLx < 45.04 


876 


8.5 X 


10- 


'5 


0.43 


44.64 


total RASS-AGN (rq) 


DR4-I- 


0.36 < 


z 


< 0.50 


44.24 < logLx < 44.99 


722 


8.3 X 


10- 


-5 


0.43 


44.61 






Optically-Selected AGNs - SDSS AGNs 














total SDSS-AGN 


DR7 


0.07 < 


z 


< 0.16 


-23.26 <Mi< -22.06 


177 








0.13 


-22.52 


total SDSS-AGN (rq) 


DR7 


0.07 < 


z 


< 0.16 


-23.01 < Mi< -22.06 


96 








0.13 


-22.45 


total SDSS-AGN 


DR7 


0.16 < 


z 


< 0.36 


-23.27 <Mi< -22.07 


3500 








0.28 


-22.55 


total SDSS-AGN (rq) 


DR7 


0.16 < 


z 


< 0.36 


-23.17 <Mi< -22.06 


2879 








0.29 


-22.51 


total SDSS-AGN (noX) 


DR7 


0.16 < 


z 


< 0.36 


-23.09 <Mi< -22.06 


2367 








0.29 


-22.47 


total SDSS-AGN (rq-(-noX) 


DR7 


0.16 < 


z 


< 0.36 


-22.98 <Mi< -22.05 


1958 








0.29 


-22.44 


total SDSS-AGN (onlyX) 


DR7 


0.16 < 


z 


< 0.36 


-23.52 <Mi< -22.11 


1133 








0.27 


-22.72 


low Mi SDSS-AGN 


DR7 


0.16 < 


z 


< 0.36 


-22.4 < Mi< -22.03 


1757 








0.28 


-22.18 


high Mi SDSS-AGN 


DR7 


0.16 < 


z 


< 0.36 


-23.60 <Mi< -22.4 


1743 








0.29 


-22.93 


low Mi SDSS-AGN (rq) 


DR7 


0.16 < 


z 


< 0.36 


-22.4 < Mi< -22.03 


1520 








0.28 


-22.18 


high Mi SDSS-AGN (rq) 


DR7 


0.16 < 


z 


< 0.36 


-23.52 <Mi< -22.4 


1359 








0.29 


-22.88 


total SDSS-AGN 


DR7 


0.36 < 


z 


< 0.50 


-23.89 < Mi < -22.36 


4404 








0.43 


-23.04 


total SDSS-AGN (rq) 


DR7 


0.36 < 


z 


< 0.50 


-23.79 <Mi< -22.36 


3773 








0.43 


-23.01 


total SDSS-AGN (noX) 


DR7 


0.36 < 


z 


< 0.50 


-23.72 <Mi< -22.35 


3421 








0.43 


-22.98 


total SDSS-AGN (rq-fnoX) 


DR7 


0.36 < 


z 


< 0.50 


-23.65 <Mi< -22.35 


2960 








0.43 


-22.96 


low Mi SDSS-AGN 


DR7 


0.36 < 


z 


< 0.50 


-22.9 < Mi< -22.22 


2059 








0.42 


-22.55 


high Mi SDSS-AGN 


DR7 


0.36 < 


z 


< 0.50 


-24.24 <Mi< -22.9 


2345 








0.44 


-23.48 


low Mi SDSS-AGN (rq) 


DR7 


0.36 < 


z 


< 0.50 


-22.9 < Mi< -22.24 


1804 








0.42 


-22.56 


high Mi SDSS-AGN (rq) 


DR7 


0.36 < 


z 


< 0.50 


-24.13 <Mi< -22.9 


1969 








0.44 


-23.43 



Note. — The notation () characterizes the average (mean) value of the given quantity. The columns 'M, Lx 
range' and '{M, Lx)' specify absolute optical magnitudes for optical samples and galactic-absorption corrected 
0.1-2.4 keV luminosities for X-ray selected samples. The absolute magnitudes are given for the SDSS filter band 
stated in the lower index and if-corrected to the redshift given in the upper index. The symbol "<" is used to 
characterize the "soft" luminosity boundary of the samples; it indicates the 10th (for the lower bound) or the 90th 
(for the upper bound) percentile. Note that unlike the main galaxy sample and the LRG sample, the extended 
LRG sample is not volume-limited. The listed co-moving number densities for the extended LRG samples are the 
number densities for the redshift range 0.36 < z < 0.42 and -23.2 < M°-^ < -21.7, where the sample is volume- 
limited. The number densities of the X-ray selected AGN samples are calculated by computing the co-moving 
volume available to each object. For optically-selected AGNs, the selection function is complex enough that we 
do not derive number densities. Abbreviations: rq - radio-quiet AGNs only (see our definition of radio-quiet in 
Sections 13.1.11 and I3.2.1|) : noX - optically-selected SDSS AGNs that are not detected by RASS; onlyX - only 
optically-selected SDSS-AGNs that are also RASS detections. 
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of M°-^ = -22.2 mag, which results in only 6292 ob- 
jects. As such a relatively low number of objects would 
yield a measured ACF with very large uncertainties, we 
use a non-volume-limited sample with —23.2 < M^-^ < 
—21.7. This is a compromise between making the sample 
volume-limited and retaining accuracy when computing 
the ACF and CCFs (see Table [I]). We plot the absolute 
magnitude versus redshift for the 28781 objects (DR7) 
in the extended LRG sample in Fig. [11 

The co-moving number density for this sample given 
in Table [1] is computed over the redshift range 0.36 < 
z < 0.42 and magnitude range -23.2 < M°-^ < -21.7, 
where the sample is volume-limited. Furthermore, it 
assumes that there is no number density evolution at 
higher redshifts (0.36 < z < 0.42) and in the range 
—23.2 < Mg-^ < —21.7. In principle, we can derive 
the co-moving number density by integrating the LRG 
luminosity function. However, the different selection 
functions for low and high redshift LRGs (cut I and II, 
lEisenstein et al.ll2d0l[ ) are difficult to model and could re- 
sult in large uncertainties. Therefore, we decide to limit 
the estimate to the volume-limited redshift range. 

The random catalogs contain 1000 times as many ob- 
jects as the data sample. The procedure follows the de- 
scription in Section [2TT2] 

3. AGN SAMPLES 

In this paper we derive the AGN ACF for X-ray and 
optically-selected SDSS AGN samples. The selection of 
the different AGN samples is described below. 

3.1. RASS/SDSS AGN Samples 

The iJOS-^T All-Sky Survey fRASS. IVoges et al.l[T999t) 
is currently the most sensitive all-sky X-ray survey. 
Although relatively shallow, this data set has a huge 
potential for science, especially when combined with 
other large-area spectrosco pic surveys such as SDSS. 
[Anderson et al.l ()2003l . [20071 ) positionally cross-correlated 
RASS sources with SDSS spectroscopic objects and 
classified RASS and SDSS-detected AGNs based on 
SDSS DR5. They find 6224 AGNs with broad permit- 
ted emission lines in excess of of 1000 kms~^ FWHM 
and 515 narrow permitted emission line AGNs match- 
ing RASS sources within 1 arcmin. More details on 
the sample sel ection are g i ven in paper I Section 2.2 
and [Anderson iEin (|2003l . [2001 . The RASS/SDSS 
AGN sample is biased toward unabsorbed AGNs due to 
ROSATs soft energy band (0.1-2.4 keV). AGNs with 
broad emission lines and UV excess are, in general, 
known to show little to no X-ray absorption and there- 
fore account for the vast majority of the RASS/SDSS 
AGNs. 

The latest available ve rsion of the RASS/SD SS AGN 
sample is based on DR5 ([Anderson et al.|[2007[ ). There- 
fore, we reconfigure the sample to publically-available 
DR4-I- geometry. We select broad emission line AGNs in 
different redshift ranges according to the redshift ranges 
of the available galaxy tracer sets (see Fig. [5] and Ta- 
ble [J). 

To study the X-ray luminosity dependence of the clus- 
tering, we further subdivide the AGNs in the redshift 
ranges of 0.07 < z < 0.16 and 0.16 < z < 0.36 
into relatively lower and higher X-ray luminosity sam- 




0.1 1.0 

redshift 



Fig. 2. — 0.1-2.4 keV observed X-ray luminosity versus redshift 
for the bro ad emiss i on hn e AGN sample in SDSS DR5 determined 
bv Anders on et al.l l|2007l '). Black symbols show objects used here 
in the redshift range of 0.07 < z < 0.50. The shaded areas show the 
different AGN subsamples, labeled with the corresponding tracer 
set for the cross-correlation measurement. The horizontal white 
lines show the dividing line between the lower X-ray luminosity and 
higher X-ray luminosity RASS/SDSS AGN subsamples for 0.07 < 
z < 0.16 log (io.l-2.4kev/[ergs-i]) = 43.54 and 0.16 < 2 < 0.36 
log (Lo.l-2.4kov/[ergs-l]) = 44.29. 

pies. The RASS/SDSS AGNs in the redshift range of 
0.36 < z < 0.50 are not subdivided into a lower and 
higher X-ray luminosity sample, as the corresponding 
tracer set (extended LRGs) has a much lower number of 
objects than the other tracer sets. Consequently, there 
are too few pairs in the CCF and prohibitively large un- 
certainties to usefully compare the RASS/SDSS AGN 
ACF of lower and higher X-ray luminosity samples. 

We use the 0.1-2.4 keV obse rved luminosity reported 
bv lAnderson et all ([2003ll2007f) . which assumes a photon 
index of F = 2.5 and is corrected for Galactic absorption. 
Our 0.1-2.4 keV luminosity cuts are log (Lx/[ergs~^]) = 
43.54 for 0.07 <z < 0.16 and log (Lx/[ergs-i]) = 44.29 
for 0.16 < z < 0.36 (see Fig. ^. Using F = 2.5, the 
luminosity cuts correspond to 0.5-10 keV luminosities of 
(Lx/[10'*'^ ergs~^]) = 1.5 and 8.5, respectively. 

A significant overlap between broad-line RASS/SDSS 
AGNs and the SDSS main galaxy sample tracer set ex- 
ists only in the lowest redshift range (0.07 < z < 0.16). 
For the low luminosity broad-line AGNs in that redshift 
range, the AGN light does not dominate the optical spec- 
trum. Hence, 186 objects of the total RASS/SDSS AGN 
sample of 629 objects are also classified as SDSS main 
galaxies (52 within the low Lx AGN sample, 134 within 
the high Lx AGN sample). In the higher redshift sam- 
ples, no match between RASS/SDSS AGN and the LRG 
tracer sets is found as the selection methods for LRGs 
and RASS/SDSS AGN are very different and high lumi- 
nosity AGNs dominate the optical spectrum. To quantify 
the effect of the overlap between both samples, we cre- 
ate a SDSS main galaxy sample that does not include 
objects that have also been classified as RASS/SDSS 
AGNs. Compared with the original CCF, the CCF of 
this new sample shows differences of less than 1% in the 
pair counts on scales larger than 1.4 Mpc. This is 
well below the statistical uncertainty of the CCF itself, 
therefore, the overlap between the AGN sample and the 
tracer set will not significantly affect our clustering mea- 
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surements in the lowest redshift range. 

We calculate the co-moving number densities given in 
Table [1] as described in detail in paper I. For a given 
R.A. and Dec we compute the limiting observable RASS 
count rate and infer the absorption-corrected flux limit 
versus survey area for RASS/SDSS AGNs. We then 
compute the co-moving volume availabl e to each object 
(Va) for being included in the sample ()Avni &: Bahcalll 
[l98m . The co-moving number density follows by com- 
puting the sum of the available volume over each object 

i 

3.1.1. Radio-quiet RASS/SDSS AGN Samples 

Radio-loud AGNs are known to be more clustered 
than radio-qu iet AGNs and res i de in very massive 
DMHs (e.g.. IMagliocc hctti ct aT 2004; Hicko x et all 
[2009t IMandelbaunTeral. 2009 ). Radio-loud AGNs are 
also more luminous in the X-rays than rad io-quiet AGNs 
(IWilkes et aLlHOOl iKnrnipe et al.ll2010bn . 

In paper I we find an X-ray luminosity dependence 
in the AGN clustering amplitude. One possible ex- 
planation is that the high ix sample contains more 
radio-loud AGNs than the low Lx sample, and there- 
fore the relative overabundance of radio-loud AGNs in 
the high Lx sample is causing the increase of the clus- 
tering amplitude. To test this hypothes is, we construct 
radio-quiet RASS/SDSS AGN samples. lAnderson et al.l 
(j200ll2007[ ) list in their table of broad-line RASS/SDSS 
AGN if an object is also detected as a radio source. 
The radio information is taken from the Faint Im- 
ages of the Radio Sky at Twen t y centimeters (FIR ST; 
IBecker. White, fc Helfandl [TOOl Iwhite et al.l Il997l) us- 
ing the NRAO Very Large Array. We therefore create 
new AGN subsamples by restricting all samples to the 
area covered by FIRST and excluding all FIRST detected 
sources, and refer to these subsamples as radio-quiet X- 
ray selected AGN samples. FIRST and SDSS cover al- 
most the same area (FIRST has ^7% less area than SDSS 
DR4+ and -5% less area than SDSS DR7). We also limit 
the tracer sets to the FIRST geometry when computing 
CCFs of radio-quiet AGN samples and the corresponding 
tracer set ACFs. 

In Tables [1] [21 and [3] we label these subsamples with 
the entry '(rq)'. Note that our approach is conservative, 
in that we do not apply the us ual radio-to-o ptical flux 
density criteria of i? > 10 (Kcll ermann et all 1989). In- 
stead we remove all radio-detected AGNs, which removes 
more objects than just those that are technically defined 
as radio-loud. Our definition of radio-quiet is that the 
AGNs are not detected by FIRST. However, for our goal 
of removing all radio-loud AGNs from the samples, the 
chosen procedure is adequate and the loss in a few ad- 
ditional AGNs will not significantly affect the clustering 
measurements and their uncertainties. 

In the redshift range of 0.07 < z < 0.16, 127 out of 504 
objects are classified as radio-quiet RASS/SDSS AGNs 
and SDSS main galaxies. Since we find no significant 
difference in the CCF of samples that include or exclude 
these objects, we use the full sample. No overlap is found 
between radio-quiet RASS/SDSS AGNs and tracer set 
objects at higher redshift s. 

In Fig. [3] we show the distribution of the RASS/SDSS 
AGNs which are flagged as radio sources by 
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Fig. 3.— Location of FIRST radio-detected RASS/SDSS AGNs 
(crosses) in the 0.1-2.4 keV observed X-ray luminosity versus red- 
shift plane for the studied DR4-I- X-ray samples. Gray dots show 
RASS/SDSS AGNs without a FIRST radio detection, which we 
refer to as 'radio-quiet' RASS/SDSS AGN. We show only those 
objects that fall within the region covered by FIRST and SDSS. 
The solid vertical lines indicate the redshift ranges of the differ- 
ent samples, while the dotted horizontal lines show the luminosity 
cuts used to create the lower and higher X-ray luminosity samples 
(see Fig. [2] for details). The percentage values give the fraction of 
radio-detected RASS/SDSS AGNs in the corresponding lower and 
higher Lx samples. 

lAnderson et~all (|2003l l2007| ). which we remove to 
create the 'radio-quiet' samples. Considering only 
objects that fall in regions covered by SDSS and FIRST, 
17% of aU RASS/SDSS AGN have radio detections in 
the 0.07 < z < 0.16 range, 10% at 0.16 < z < 0.36, 
and 14% at 0.36 < z < 0.50. We also report the 
radio-detected AGN fraction in the corresponding lower 
and higher Lx samples as the percentage values listed 
in Fig. H 

3.1.2. Narrow-line RASS/SDSS AGN Samples 

In order to test narrow-line versus broad-line AGN 
clustering, we construct X-ray selected narrow-line AGN 
samples . In addition t o the predominant broad-line 
AGNs, lAnderson etall (l2007h classified -7% of aU 
RASS/SDSS AGNs as X-ray emitting AGNs having nar- 
rower permitted emission lines. These 515 objects consist 
of X-ray emitting AGN subclasses such as narrow-line 
Seyfert I's (NLSls) (10%), Seyfert 1.5 (29%), Seyfert 1.8 
(18%), Seyfert 1.9 galaxies (22%), and Seyfert 2 candi- 
dates (21%). In total, 22% of these objects also have a 
FIRST radio detection. NLSls, Seyfert 1.5, 1.8, and 1.9 
AGNs show a mix of narrow and broad permitted typ- 
ical AGN line components, while Seyfert 2 candidates 
have only narrow permitted emission lines. The latter 
are called candidates because it is well known that a 
fair fraction of them turn out to be reclassified as NLSl, 
Seyfert 1.8, or Seyfert 1.9 when re- observed with signifi - 
cantly improved spectroscopy (e.g.. lHalpern et a"I]|1999D . 
Seyfert 1.5 and 1.8 galaxies have optical spectra with a 
broad line II/3 component (exceeding FWHM values of 
2500 kms~^) at a very low fiux level. For more details 
about the different narrow - line AGN s ubclasses in this 
sample see lAnderson et all ()2003l . l2007D . Figure S] shows 
that the narrow- line RASS/SDSS AGNs are mainly iden- 
tified only at lower redshifts, to z ^ 0.35. Therefore, 
their classification as narrow-line AGNs is based on the 



Main Galaxy 
sample 



21^ 



12% 



15% 



XX„KXX X,? 

, ..X x^^/'^x'^x 

^111 ?fl ^111 lX^7l IXIXI ^ I I I M 

X Sj,^ x< / xX 
x'' XX X 

• v ■^>^. ~ 

■' 7% 



XX X 
, XX X 



LRG sample 



T3 P 



8 



Krumpe et al. 



permitted H/3 and Ha lines. 

The studied narrow-line RASS/SDSS AGNs are found 
to have, on average, lower observed X-ray luminosities 
than broad-line RASS/SDSS AGNs. However, the vast 
majority of broad-line RASS/SDSS AGNs are known to 
be unabsorbed X-ray sources, therefore, their observed 
X-ray luminosity is equal to their intrinsic X-ray luminos- 
ity. On the other hand, narrow-line RASS/SDSS AGNs 
may be absorbed AGNs and their intrinsic X-ray lumi- 
nosities could be higher than the observed ones given 
in Figure m Consequently, both classes of RASS/SDSS 
AGN samples may be very similar with respect to intrin- 
sic X-ray luminosity. 

Due to ROSATs soft energy range of 0.1-2.4 keV, 
narrow-line RASS/SDSS AGNs can be absorbed only 
by moderate column densities (Nn). We simulate an 
X-ray spectrum with a photon index of F ~ 2.5, galac- 
tic absorption of N^^gai = 2 x 10^°cm~^, and a typ- 
ical redshift of z = 0.15, in order to estimate which 
intrinsic A^h values are detectable with ROSAT. Com- 
pared to an unabsorbed source, the flux in the 0.1-2.4 
keV ROSAT band drops down to 61% when an intrin- 
sic column density of A'h = lO'^^ cm~^ is used (21% 
for A'h = lO^^cm"^). We conclude that only objects 
with A^H < 10^2 cm-2 are detected in the RASS. Most 
likely the narrow-line RASS/SDSS AGN sample consists 
mainly of a mixture of unabsorbed and only moderately 
absorbed (a few A^h = 10^^ cm'^) AGNs. 

Although broad-line AGNs have a much lower pro- 
jected space density than narrow-line AGNs and there- 
fore yield, in general, more reliable identification with 
the RASS counterpart. I Anderson et al.l ()2003ll2007 D suc- 
cessfully demonstrate that the narrow-line RASS/SDSS 
AGN are statistically very reliable identifications as well. 
They estimate that less than 5% of the counterparts are 
spurious random chance positional coincidences. Con- 
vincing evidence derives from the distribution of posi- 
tional offsets relative to the X-ray positional error and 
equal area annuli. Furthermore, the observed distribu- 
tion of the ratios of RASS/SDSS X-ray-to-optical flux 
matches expect ations for typica l X-ra y emit ting AGNs. 
For details see [Anderson etall (|2003l 120071 ). We con- 
clude that the narrow-hne RASS/SDSS AGNs have a 
comparable high counterpart reliability as the broad-line 
RASS/SDSS AGNs. 

In the redshift range of 0.07 < z < 0.16 31 out of 
194 narrow-line RASS/SDSS AGNs are also classified 
as SDSS main galaxies, while at 0.16 < z < 0.36 only 
one object belongs to the narrow-line RASS/SDSS AGN 
sample and the LRG tracer set. As for the broad-line 
RASS/SDSS AGNs, we create a SDSS main galaxy sam- 
ple that does not include galaxies that have been also 
classified as narrow- line RASS/SDSS AGNs and com- 
pute the corresponding CCF. On scales greater than 5 
Mpc, the pair counts differ by less than 0.1% from 
the original CCF. Therefore, the overlap in the samples 
does not affect the clustering results for the narrow-line 
RASS/SDSS AGNs. 

3.2. Optically-selected SDSS Samples 

All of our optically-selected SDSS AGNs (called 
'quasars' in the SD SS literature) are drawn from 
[Schneider erall (|2010l ) and use the full SDSS survey 
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Fig. 4.— Location of narrow-line RASS/SDSS AGNs (black 
squares) in the 0.1-2.4 keV observed X-ray luminosity versus red- 
shift plane for the studied DR4-(- X-ray samples. Gray dots show 
broad-line RASS/SDSS AGNs. The filled areas illustrate the red- 
shift ranges of the different samples used here. 

(DR7) . The AGN candidate selection is described in de- 
tail in [Richards et al] (|2002D and can be summarized as 
follows. The highest priority is given to FIRST-detected 
optical point sources. Then sources with non-stellar col- 
ors in the ugriz photometry data are considered. Objects 
with photometric redshifts of z < 3 are targeted even 
if they are spatially resolved. FIRST-detected optical 
point sources and z < 3 AGN candidates are required to 
have a Galactic extinction-corrected i magnitude of 19.1. 
This selection method picks ~ 18 objects per square de- 
gree which are followed up with the SDSS spectrograph. 
The resulting 'primary' SDSS AGN sample is then sup- 
plemented by objects targeted by other SDSS spectro- 
scopic selections (main galaxies, LRGs, RASS, stars, and 
serendipitous sources) that turned out to be AGNs. In 
the redshift range of interest (0.07 < z < 0.50), these sec- 
ondary channels account for ^9% of the total SDSS AGN 
sample (vast majority from the galaxy target selection). 
If a significant number of SDSS AGNs from these sec- 
ondary channels were originally selected from the LRG 
sample, which is known to be strongly clustered, it could 
bias the AGN clustering. We verified that there is no 
overlap between the SDSS AGN and the LRG sample. 
Therefore, we do not expect any infiuences on our AGN 
clustering measurements caused by the secondary selec- 
ti on channel. 

[Schneider et al.l (|2010f) (and references within) have 
constructed AGN catalogs based on different SDSS data 
releases. They apply a luminosity selection of Mi < — 22 
mag and require that objects have at least one emis- 
sion line exceeding a FWHM of 1000 kms~^. Ob- 
jects that have a spectrum with only narrow permitted 
AGN-typical emission lines are removed. The absolute 
magnitude Mi is computed by using the j PSF Galac- 
tic extinction-corrected magnitude measurement and as- 
suming a typical AGN spectral energy slope. 

Historically, quasars are defined as objects at the high 
end of the AGN luminosity functi on having Mb < —23 
mag fe.g.. [Schmidt fc Green[[T983h . M, = -22 mag cor- 
responds to Mb = —22.4 mag for an typical AGN at 
z = 0. The Schneider et al. SDSS AGN catalog papers 
use the i-band instead of the i3-band in part because the 
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z-filter is less affected by Galactic absorption. However, 
a significant disadvantage of this is that the host galaxy 
light may represent a larger contribution to the total flux 
then the AGN. AGNs near the Mi = -22 mag cut can 
be equally bright as the host galaxy, e.g., host galaxies 
at 2: ~ 0.4 with i = 19.1 mag (within the detection limit 
of the SDSS AGN selection method) have Mi = -22 
mag. Consequently, these AGNs may be less luminous 
than their quoted optical magnitude. This effect is some- 
what mitigated by the use of the PSF photometric data 
in the AGN selection but should be kept in mind when 
interpreting our results. 

The AGN candidates selection by iRichards et al.l 
(|2002f ) has und ergone constant modif i cation to improve 
the efficiency. [Schneider eFall (|20ia 120071 ) use the fi- 
nal selection algorithm. AGNs in these catalogs have 
two spectroscopic target selection flags: BEST (final al- 
gorithm) and TARGET (actual algorithm during target- 
ing). BEST uses the latest photometric software and has 
the highest quality data. The continuous modification of 
the AGN selection method and the inclusion of AGNs 
not selected by the standard selection means that the 
AGN cat alogs are not sta t istica lly clean samples. The 
catalog of iSchneider et al.1 ((2010!) contains 105,783 spec- 
troscopically confirmed AGNs which have luminosities of 
Mi < —22, have at least one emission line exceeding a 
FWHM of 1000 kms~^, have highly reliable redshifts, 
and are fainter than i ^ 15. We extract only objects 
that have a MODE flag of 'PRIMARY'. This procedure 
applies to 99% of all objects and limits our sample to ob- 
jects that have been spectroscopically followed up based 
on a target selection and are not blended. We use the 
BEST flags whenever TARGET and BEST are available. 
We use optically-selected SDSS AGNs in the redshift 
ranges of the corresponding tracer sets (0.07 < z < 0.16, 
0.16 < z < 0.36, 0.36 < z < 0.50). Except for the 
0.07 < z < 0.16 AGN sample, where the number of ob- 
jects is very low, we further subdivide the samples into 
lower and higher Mi subsamples (see Table [Ij. The Mi 
cuts are chosen to yield approximately the same number 
of objects in the different luminosity subsamples. We 
use Mi = -22.4 for the 0.16 < z < 0.36 sample and 
M, = -22.9 for the 0.36 <z < 0.50 sample (Fig.[S]). 

In the 0.07 < z < 0.16 redshift rang e we find a high 
overla p between objects derived from the lSchneider et all 
(poTol ) AGN sample and the SDSS main galaxy sample. 
133 out of 177 SDSS AGNs are also classified as SDSS 
main galaxies. As with the RASS/SDSS AGNs, we ex- 
clude these objects from the SDSS main galaxy sample 
and compute the CCF. The difference in the CCF oc- 
curs only on the larger scales and is less than 1%; there- 
fore the overlap affects the clustering measurements very 
little. At higher redshifts, there is no overlap between 
optically-selected SDSS AGN samples and SDSS LRGs 
used as our tracer sets. For the optically-selected SDSS 
AGNs we do not estimate the co-moving number densi- 
ties in Table [TJ Their selection function is very complex 
and requires a sophisticated modeling method, which is 
beyond the scope of this paper. 

3.2.1. Radio-qmet Optically-selected SDSS AGN Samples 

To test the impact of radio-loud AGNs on the clus- 
tering signal of optically-selected AGN samples, we 
further exclude FIRST radio-detected objects in the 




II 



iir 
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Fig. 5. — Absolute magnitude Mi vs. redshift for optically- 
selected broad-line SDSS DR7 AGNs from [Schneider et al. (2010). 
The vertical black lines illustrate the redshift ranges of the different 
AGN samples and tracer sets, foUowring Figs. [2] and [4] The white 
horizontal solid lines show the cut used to create lower and higher 
Mi optically-selected SDSS AGN subsamples. 




0.1 0.2 0.3 0.4 0.5 

redshift 

Fig. 6. — Similar to Figure \5\ here showing the location of 
optically-selected broad-line FIRST-detected AGNs (top panel, 
crosses) and RASS-detected AGNs (lower panel, asterisks) in the 
absolute Mi-magnitude-redshift plane. In the top panel we select 
only objects that fall in regions covered by FIRST and SDSS. The 
solid vertical lines show the different redshift ranges of the samples. 



ISchneider et al.l ()2010D AGN samples. An AGN has a 
radio detection if its position coincidences with a FIRST 
catalog entry within 2 arcsec. We restrict all samples 
to the FIRST area to ensure that no radio-loud AGN is 
selected for the radio-quiet samples. Then we remove ob- 
jects that have a measured FIRST radio flux in the AGN 
catalog. This is the same conservative approach of ex- 
cluding aU FIRST-detected AGNs as we do in Sect. l5XT] 
for the radio-quiet RASS/SDSS AGN subsamples. Here 
we also use our deflnition of radio-quiet AGN samples 
for the optically selected AGNs, namely that radio-to- 
optical flux density R based on the FIRST is i? = 0. We 
also account for the restriction to the FIRST geometry 
in the corresponding tracer sets when we derive the ACF 
of the radio-quiet optically-selected AGN samples. The 
resulting samples are labeled with the subsequent entry 
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'(rq)' in Tables andUl 

Figure |6] (top panel) shows the distribution of 
optically-selected FIRST-detected AGNs within the 
larger optically-selected AGN sample. A fraction of 44% 
of the 0.07 < z < 0.16 AGNs, 13% of 0.16 < z < 0.36 
AGNs, and 9% of the 0.36 < z < 0.50 AGNs are detected 
by FIRST (considering only regions that are covered by 
FIRST and SDSS). The fraction is a sensitive function 
of redshift, as it depends on luminosity. The higher 
subsamples contain more FIRST radio-detections (17% 
in 0.16 < z < 0.36 and 12% in 0.36 < z < 0.50) than 
the lower Mi subsamples (8% in 0.16 < z < 0.36 and 
5% in 0.36 < z < 0.50). 78 out of the 96 radio-quiet 
SDSS AGN in the redshift range 0.07 < z < 0.16 are 
also classified as SDSS main galaxies. 

3.2.2. X-ray selected Optical AGN Samples 

The large number of X-ray detected optically-selected 
SDSS AGNs at redshifts above z = 0.16 allows us to 
test how different selections of AGNs affect the cluster- 
ing signal of broad-lin e AGNs. The catalog paper by 
[Schneider et al.l (j201 0D lists for each individual AGN the 
relevant information if the object has a detection in the 
RASS faint or bright source catalog. Beside the con- 
struction of a radio-quiet optically-selected SDSS AGN 
sample (Section I3.2.ip . we create a sample that con- 
tains only optically-selected broad-line SDSS AGNs that 
are not detected in RASS. This sample has the subse- 
quent entry '(noX)' in the corresponding tables. FigurelH] 
(panel b) shows the distribution of objects with RASS 
detections among the optically-selected SDSS AGN sam- 
ple. A fraction of 70% of the optically-selected SDSS 
AGNs are also detected by RASS in the redshift range 
of 0.07 < z < 0.16. At higher redshifts RASS detects 
32% (0.16 <z< 0.36) and 22% (0.36 <z < 0.50) of the 
optically-selected broad-line SDSS AGNs. This is not 
surprising, as the low ROSAT sensitivity results in the 
RASS redshift distribution rising quickly to z ~ 0.15 and 
strongly decreasing at higher redshift. 

The optical and X-ray luminosities of broad emission 
line AGNs are known to be strongly correlated, the ratio 
of which is often expressed by the optical-to-X-ray spec- 
tral energy i ndex apx (e.g.. lAvni fc TananbaumI 119861: 



Green et al.1 [19951: ISteffen et al.l 120061: iKrumpe et alj 



2007t [Anderson et al. 112003). Therefore, RASS detects 



44% of the AGNs in the higher Mi sample at 0.16 < 
z < 0.36 (29% at 0.36 < z < 0.50) but only 24% of the 
AGNs in the lower Mi sample at 0.16 < z < 0.36 (17% 
at 0.36 < z< 0.50). 

RASS contains only the highest flux X-ray emitting 
AGNs but has the complementary advantage of de- 
tecting lower AGN activity compared to the optically- 
selected broad-line SDSS AGNs. X-ray luminosities of 
log (Lx/[ergs~^]) > 42 are a clear indicator of AGN 
activity, while in the optical a strong starlight compo- 
nent from the host galaxy can make it difficult to detect 
broad-line AGNs. The location of the RASS-detected 
AGNs in Fig. [6] (panel b) indicates that RASS extends 
the detections of broad-line AGNs below the o ptical cut 
of M, = -22 mag used by [Schneider et al.l (|2010f ) at 
z < 0.35. 

Furthermore, we also create AGN samples that con- 
tain only optically-selected SDSS AGN that are also de- 
tected as RASS sources (Fig. [6j panel b) . The samples 



are labeled with the subsequent entry '(onlyX)'. For these 
object we have both the X-ray and optical luminosities 
{Lx, Mi). In general, as expected by the Oox connection, 
high/low Mi corresponds to high/low Lx- The fraction 
of FIRST-detected AGNs increases more with Mi than 
with Lx- 

Finally, we design samples of optically-selected SDSS 
AGNs that are neither FIRST nor RASS detections. 
These samples only cover the FIRST area. This selec- 
tion results in 56% of the 0.16 < z < 0.36 and 67% of the 
0.36 <z< 0.50 optically-selected SDSS AGNs. We refer 
to these samples by the subsequent entry '(rq-fnoX)'. 

4. CLUSTERING ANALYSIS 

We measure the two-point correlation function ^(r) 
(jFecblcs 1980), which measures the excess probability 
dP above a Poisson distribution of finding an object 
in a volume element dV at a distance r from another 
randomly chosen object. The auto-correlation function 
(ACF) measures the spatial clustering of objects in the 
same sample, while the cross-correlation function (CCF) 
measures the clustering of objects in two different sam- 
ples. We use the same approach as described in detail 
in paper I (Section 3). Here we reiterate the essential 
elements of our method. 

We use the correlation estimator of iDavis fc Peebles! 
(|1983D in the form 



DD{r) 
DR{r) 



1 



(2) 



where DD{r) is the number of data-data pairs with a sep- 
aration r, and DR{r) is the number data-random pairs; 
both pair counts have been normalized by the number 
density of data and random points. We measure ^ on a 
two-dimensional grid of separations r^, perpendicular to 
the line of sight, and tt, along the line of sight, to separate 
the effects of redshift space distortion due to peculiar ve- 
locities along the line of sight. We obtain the projected 
correlation function Wp{rp) by integrating (,{rp,Tr) along 
the TT direction. 

As in paper I, we infer the AGN ACF from the CCF 
between AGNs and correspondi ng galaxy tracer set and 
the ACF of the tracer set using iCoil et al.l (|2009[ ): 



Wp{AGN\AGN) 



[wp{AGN\TRACE)Y 
Wp{TRACE\TRACE) 



(3) 



where Wp{AGN\AGN), Wp{T RAG E\T RAG E) are the 
ACFs of the AGN and the corresponding tracer set, re- 
spectively, and Wp{AGN\TRAGE) is the CCF of the 
AGNs with the tracer set. In other words, we assume 
that the CCF is the geometric mean of two ACFs, which 
has been verified to be valid by Zcha vi~et al.l ()2011[ ) (Ap- 
pendix A). 
The CCF is computed by applying Eq. [2] 



N- TRACE 



-Pagn -Ptrace 
-Dagn -Rtrace 



1. 



(4) 



For our purposes, the use of this simple estimator has 
several major advantages and results in only a marginal 
loss in the signal-to-noise r atio when compared t o more 
advanced estimators (e.g., iLandv fc Szalavl [19931) . The 
estimator in Eq. |4] requires the generation of a random 
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catalog only for the tracer set. The tracer sets have well- 
defined selection functions and are, except for the ex- 
tended LRG sample, volume-limited. Since the random 
catalog should exactly match all observational biases to 
minimize the systematic uncertainties, well understood 
selection effects are a key to generation proper random 
samples. The AGN samples suffer from very complex and 
hard to model selection functions. Therefore, a random 
catalog of X-ray selected RASS/SDSS AGN is subject 
to large systematic uncertainties due to the difficulty in 
accurately modeling the position-dependent sensitivity 
limit and the variation in the fiux limit of the sources 
(caused by changing Galactic absorption over the sky 
and spectrum-dependent corrections) . Optically-selected 
SDSS AGNs (see Section 13. 2[) rely on constantly modi- 
fied selection algorithms and the acceptance of additional 
incomplete AGN selection methods. Consequently, the 
modeling of their selection function for the generation of 
a random catalog would be very challenging. 

The errors in the adjacent bins in correlation measure- 
ments are not independent. Poisson errors will signif- 
icantly underestimate the uncertainties and should not 
be used for error calculations. Instead, we use the jack- 
knife resampling technique to estimate the measurement 
errors as well as the covariance matrix Mij , which reflects 
the degree to which bin i is correlated with bin j. The 
covariance matrix is used to obtain reliable power law 
flts to Wp(rp) by minimizing the correlated values. 
In our jackknife resampling, we divide the survey area 
into A't = 100 subsection for the DR4-|- geometry and 
131 subsections for DR7, each of which is ~50-60 deg^. 
These jackknife-resampled correlation functions de- 
fine the covariance matrix (Eq. [5]) : 



iVx - 1 



.fc=l 



X Wfe(rpj) - {w{rp,j)) 



(5) 



We calculate Wp{rp) times, where each jackknife sam- 
ple excludes one section and w^irp^i) and Wk{rpj) are 
from the k-th jackknife samples of the AGN ACF and 
{w{rp^i)), {w{rpj)) are the averages over all of the jack- 
knife samples. The uncertainties represent a la (68.3%) 
confidence interval. 

The generation of covariance matrix for the inferred 
AGN ACF considers the A^t jackknife-resampled corre- 
lation functions of the CCF (AGN and corresponding 
tracer set) and the tracer set ACF. For each jacknife sam- 
ple we calculate the inferred AGN ACF by using Equa- 
tionO The resulting Wp(rp) jackknife-resampled pro- 
jected correlation functions of the inferred ACFs are then 
used to compute the covariance matrix of the inferred 
AGN ACF. 

4.1. Inferring the AGN Auto-correlation Function 

To infer the AGN ACF, we measure the CCF of the 
AGN sample with the tracer set and the ACF of the 
tracer set. In both cases we measure rp in a range of 
0.05-40 Mpc in 15 bins in a logarithmic scale. The 
upper 11 bins are identical with the bins used in paper I 
and cover the rp range of 0.3-40 Mpc. Consequently, 
we extend the measurements by four additional bins to 
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Fig. 7. — Contour plots of the auto-correlation functions ^{vp, it) 
of the different tracer sets (left panels, top - SDSS main galaxies, 
middle - LRGs, bottom - extended LRGs) and cross-correlation 
functions of the X-ray selected (RASS/SDSS) AGN samples (mid- 
dle panels) and the optically-selected SDSS AGNs (right panels) 
with the corresponding tracer set. The shown CCFs use in all cases 
the total AGN sample in the corresponding redshift range. In the 
TT direction we use a binning of 5 Mpc. Contour lines show 
constant correlation strength for the two-dimensional correlation 
function ^{rp, tt). The contour levels are 0.0 (dotted line), 0.1, 0.2, 
0.5, 1.0 (thick solid line), 2.0, and 5.0. 

smaller scales. We compute tt in steps of 5 Mpc in 
a range of tt = — 200 Mpc. The resulting £,{rp, tt) 
are shown in Fig.[7]for the ACF of the tracer sets and for 
the CCF of the total AGN samples (X-ray & optically- 
selected) with the corresponding tracer sets. 

Although the projected correlation function is com- 
puted by integrating over tt to infinity (see paper I, Equa- 
tions 5, 11), in practice an upper bound of the integration 
(TTmax) IS used to includc most of the correlated pairs, 
give stable solutions, and suppress the noise introduced 
by distant, uncorrelated pairs. We compute Wp{rp) for a 
set of TTmax ranging from 10-160 Mpc in steps of 10 
ft-^^ Mpc. We then fit Wp{rp) over a rp range of 0.3-40 
h^^ Mpc with a fixed 7 — 1.9 and determine the correla- 
tion length ro for the individual TTmax measurements. As 
in paper I, we find that the LRG ACFs (LRG sample and 
extended LRG sample) saturate at TTmax = 80 Mpc. 
All CCFs and the main galaxy sample ACFs saturate 
at TTmax — 40 Mpc. In addition, above these val- 
ues the corresponding correlation lengths do not change 
by more than Icr, considering the increased uncertain- 
ties with increasing TTmax values. The use of TTmax = 80 
h-^ Mpc for the LRG ACFs and TTmax = iO h'^ Mpc 
for the main galaxy ACFs and all CCFs matches the 
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Ti'max values used for these samples by other stud i es, e.g. , 
IZehavi et al.l (l2005aD (LRG) and IZehavi et al.l (I2005b[ ) 
(SDSS galaxies). The Wp{rp) CCFs for the different to- 
tal AGN samples with the corresponding tracer sets are 
shown in Fig.[51 while the resulting power law fits for the 
ACFs and CCFs based on 



Wp{rp 



roY r(l/2)r((7^1)/2) 

r(7/2) 



(6) 



where r(a;) is the Gamma function, are listed in Table[2l 
Our values for the ACFs of the tracer sets agree well 
with measurements from other studies. Using a slightly 
different magnitude cut of —22 < Mr < —2 1 in the 
redshift range 0.07 < z < 0.16, IZehavi et all (|2005bf ) 
find for the SDSS DR2 main galaxy sample a correla- 
tion length of ro = 6.16 ± 0.17 Mpc and a slope 
of 7 = 1.85 ± 0.03, fitting over the range rp=0. 13-20 
h^^ Mpc. IZehavi et al.l (|2005al ) study the clustering of 
~30,000 LRGs with -23.2 < M°-^ < -21.2 and measure 

ro = 9.80 ± 0.20 h'^ Mpc and a slope of 7 = 1.94 ± 0.02, 
fitting over the range rp=0.3-30 h"^ Mpc. Their cluster- 
ing measurement in a sample of high luminosity LRGs 
(-23.2 < M°-3 < -21.8, 0.16 < z < 0.44) yields 

ro = 11.21 ± 0.24 /i-i Mpc and 7 = 1.92 ± 0.03. 

In Table |3] we list the redshift range, the median ef- 
fective redshift of Nccf{z) for the corresponding AGN 
samples, the derived best ro and 7 values based on power 
law fits, and ro for a power law fit with a fixed slope of 
7 = 1.9. The data are fitted over the range rp = 0.3 — 15 
h^^ Mpc to be consistent with paper I. Since we measure 
the CCF to infer the ACF, the resulting effective red- 
shift distribution for the clustering signal is determined 
by both the redshift distribution of the tracer set and the 
AGN sample: A^ccf(^) = A^tracer(2:) * Nagn{z). 

The clustering strength is commonly expressed in 
terms of the rms fluctuations within a sphere with a co- 
moving radius of 8 Mpc (cr8,AGN, see Equation 13 
in paper I). We derive (Ts.agn from the best fit pa- 
rameters of the power law fits. The uncertainties on 
fs.AGN are derived from the ro versus 7 confidence con- 
tours of the one-parameter fit based on a correlated 
= Xmin + 1-0- Using fT8,AGN, we further derive the bias 
parameter b = crs,AG'Niz)/as{z) based on our power law 
fits and give these values in column 'b{z) PL-fits' of Ta- 
ble [21 The parameter b indicates the clustering strength 
by comparing the observed AGN clustering to that of the 
underlying mass distribution from linear growth theory 
()Hamilton|[200l . with as{z) = D{z)asiz = 0), where 
D{z) = D^(z) I D^(z = O) is the linear growth factor 
(Section 7.5 in lDodelson|[200l . We use a^iz = 0) = 0.8 
(see Section [T|). The bias quantifies the amplification fac- 
tor of the contrast of the object distribution with respect 
to that of the dark matter density distribution. The un- 
certainties of b are derived from the standard deviation 

of crg.AGN- 

Column '6(2) HOD' of Table [3] lists the bias param- 
eter derived using halo occupation distribution (HOD) 
modeling, which is described below in Section [S] Us- 
ing Equation 8 of iSheth et al.l (I2001D and the impr oved 
fit for this equation given by iTinker et al.l (12003) . we 
compute the expected large-scale Eulerian bias factor 
for different DMH masses at different redshifts. Com- 



TABLE 2 

Power Law Fits to the ACFs of the Tracer Sets and the 
CCFs OF THE AGN - Tracer sets 



Sample 


Redshift 


ro 


7 


Name 


Range 


(/i-l Mpc) 





main galaxy (DR4-I-) 
main galaxy (DR7) 
LRG (DR4-(-) 
LRG (DR7) 
extended LRG (DR4-I-) 
extended LRG (DR7) 



SDSS Tracer Sets 
0.07 - 0.16 
0.07 - 0.16 
0.16 - 0.36 
0.16 - 0.36 
0.36 - 0.50 
0.36 - 0.50 



,+0.12 
-0.12 
,+0.12 
-0.12 

,+0.14 
-0.14 
+0.13 
-0.13 
,+0.24 
0.24 



6.30: 

6.27j 
9.63j 
9.54j 
10.90:! 



X-ray Selected AGN - RASS/SDSS AGN 



total RASS-AGN 0.07 - 0.16 

total RASS-AGN(rq) 0.07-0.16 

low Lx RASS-AGN 0.07 - 0.16 

high Lx RASS-AGN 0.07-0.16 

low Lx RASS-AGN (rq) 0.07 - 0.16 

high Lx RASS-AGN (rq) 0.07 - 0.16 

narrow line RASS-AGN 0.07-0.16 

total RASS-AGN 0.16 - 0.36 

total RASS-AGN(rq) 0.16 - 0.36 

low Lx RASS-AGN 0.16 - 0.36 

high Lx RASS-AGN 0.16-0.36 

low Lx RASS-AGN (rq) 0.16 - 0.36 

high Lx RASS-AGN (rq) 0.16 - 0.36 

narrow line RASS-AGN 0.16-0.36 

total RASS-AGN 0.36 - 0.50 

total RASS-AGN(rq) 0.36 - 0.50 



5.79 
5.78 
5.19 
6.07 
5.49 
6.06 
4.99 



+0.24 
0.25 
+0.27 
0.27 
+0.34 
0.36 
+0.34 
0.35 
+0.35 
0.36 
+0.42 
0.44 
+0.45 
0.48 



i+0.27 
-0.28 
7 m+0-28 
'■'-'J^-0.29 
-,+0.36 
-0.38 
-1+0.38 
-0.41 
-,+0.37 
-0.39 
3+0.39 
-0.41 
a+0.74 
-0.86 



6.88^ 



6.22; 
7.50j 
6.09j 
7.83j 
5.78^ 



6.82^ 



^+0.42 
-0.44 

"■"'-0.49 



Optically-selected AGN - SDSS AGN 



total SDSS-AGN 0.07 - 0.16 4.93 

total SDSS-AGN(rq) 0.07 - 0.16 4.70 

total SDSS-AGN 0.16 - 0.36 6.91 

total SDSS-AGN(rq) 0.16 - 0.36 6.98 

total SDSS-AGN (noX) 0.16 - 0.36 6.93 

total SDSS-AGN (rq+noX) 0.16 - 0.36 6.96 

total SDSS-AGN (onlyX) 0.16 - 0.36 7.10 

low Mi SDSS-AGN 0.16 - 0.36 6.91 

high Mi SDSS-AGN 0.16 - 0.36 6.74 

low Mi SDSS-AGN (rq) 0.16 - 0.36 

high Mi SDSS-AGN (rq) 0.16 - 0.36 6.91 

total SDSS-AGN 0.36 - 0.50 

total SDSS-AGN (rq) 0.36 - 0.50 

total SDSS-AGN (noX) 0.36 - 0.50 

total SDSS-AGN (rq+noX) 0.36 - 0.50 

low Mi SDSS-AGN 0.36 - 0.50 7.12 

high Mi SDSS-AGN 0.36 - 0.50 6.96 

low Mi SDSS-AGN (rq) 0.36 - 0.50 7.12 

high Mi SDSS-AGN (rq) 0.36 - 0.50 7.16 



+0.41 
0.43 
+0.53 
0.58 

+0.17 
0.18 
+0.20 
0.20 
+0.21 
0.22 
+0.23 
0.23 
0.32 
0.33 
+0.24 
0.25 
0.22 
0.22 
+0.26 
-0.26 
+0.27 
0.28 

7 n-i+0.21 
' -^^^-0.22 

' -^^-0.23 

' •'-'^-0.25 

7-o8iE;:^| 

+0.30 
0.31 
+0.33 
0.36 
+0.32 
0.34 
+0.33 
0.36 



1.85 
1.84 

1.98 
1.95 
1.91 
1 



+0.02 
0.02 
+0.02 
0.02 
+0.02 
0.02 
+0.02 
0.02 

+0.03 
0.03 
+0.03 
-0.02 



1.84 
1.86 
1.96 
1.79 
1.94 
1.79 
1.82 

1.85 
1.86 
1.85 
1.90 
1.85 
1.94 
1.77 



+0.05 
0.05 
+0.05 
0.05 
+0.10 
0.09 
+0.06 
0.06 
+0.09 
0.09 
+0.06 
0.06 
+0.12 
0.10 

0.04 
0.04 
+0.05 
0.05 
+0.07 
0.06 
+0.07 
0.07 
+0.07 
0.06 
+0.08 
-0.08 
-0.14 
-0.14 



1.97^ 



fO.lO 
-0.09 



1.98 
2.11 



+0.12 
0.11 

+0.17 
0.14 



1.91'' 



(-0.04 
-0.04 



-0.04 
,+0.04 
-0.04 



1.93^ 



1.85 
1.91 
1.90 
1.90 
1.87 



1.87 
1.91 
1.81 
1.84 
1.91 
1.81 
1.95 
1.84 



0.04 
+0.05 
0.05 
+0.05 
0.05 
+0.05 
0.05 
+0.05 
0.05 
+0.05 
-0.05 

-0.04 
-0.04 
+0.04 
0.04 
+0.04 
0.04 
+0.04 
0.04 
+0.06 
0.06 
+0.07 
0.07 
+0.06 
0.06 
+0.06 
0.06 



Note. — Values of ro and 7 obtained from a power 
law fit to Wp(rp) over the range rp=0.3-40 Mpc for all 
samples using the full covariance matrix. For the LRG and 
extended LRG ACFs, we use TTmax = 80 h^^ Mpc, while 
for all other ACFs and CCFs we use TTmax = 40 Mpc. 
See Tabled] for the definition of the samples. 
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paring the observed h value from HOD modeling with 
the DMH bias factor from ACDM cosmological simu- 
lations provides an estimate of the typical DMH mass 

(&dmh(Mdmh) = &0BS30d) in which the different AGN 
samples reside, listed in the last column of Tabled 

4.2. Robustness of the Clustering Measurements 

In this section, we verify the stability of our result 
against possible observational biases and systematic ef- 
fects. As shown in paper I, the somewhat non-contiguous 
coverage of the SDSS DR4-I- survey does not influence the 
clustering results significantly, given the uncertainties. 
For DR7 the situations improves as the SDSS geometry 
is much more contiguous. We verify that the number of 
random points is large enough to lead to a high number 
of pair counts at the smallest scales measured and not 
i ntroduce noise . 

iZehavi et al.l ()2005b[ ) note that the largest structure 
detected in SDSS (the Sloan Great Wall) influences the 
galaxy clustering significantly for samples with Af"'^ < 
—21 and z < 0.1. Therefore, we explore this effect on 
our SDSS main galaxy sample. We measure the ACF 
for a main galaxy sample that excludes the Sloan Great 
Wah (165 < RA < 210 and -5 < Dec < 5). We find 
a correlation length of rg = 6.35 ± 0.12 Mpc and 
7 = 1.83 ± 0.02, which agrees well with the clustering 
measurements in which the Sloan Great Wall is included 
(see Table [2|). We conclude that our main galaxy ACFs 
and the CCFs using the main galaxy sample as a tracer 
set are not affected by this supercluster at z ^ 0.08. 

For several of the AGN subsamples split by luminosity, 
we have tested that slightly changing the luminosity cuts 
by up to ±0.2 mags (both brighter and fainter) does not 
significantly change the measured CCFs. The combina- 
tion of the different tests listed above provides convincing 
evidence that our results are not significantly infiuenced 
by systematic effects and demonstrates their robustness. 

5. BIAS FROM THE HOD MODELING 

In paper II, we develop a novel method to infer the halo 
occupation distribution (HOD) of RASS/SDSS AGNs di- 
rectly from the well-constrained cross-correlation func- 
tion of RASS/SDSS AGNs with LRGs. The resuhs from 
paper II show that the linear bias parameters and typ- 
ical DMH masses derived from the best power law fits 
down to Tp « 0.3 h^^ Mpc are subject to systematic er- 
rors. This is mainly because the power law fits include 
scales in the non-linear regime (rp < 1.5 h^^ Mpc), where 
the contribution from pairs of objects that belong to the 
same DMH (the one halo term) is substantial. In this 
non-linear regime the bias-DMH mass relation based on 
linear theory, in principle, should not be applied. An- 
other, less significant source of systematic error is that 
even in the linear regime (rp > 1.5 h^^ Mpc), the under- 
lying matter correlation function deviates from a power 
law. In order to avoid these issues, lAUevato et all poTTI ) 
derive the bias parameters of their AGN samples in the 
XMM-COSMOS survey by modeling their AGN ACFs 
at rp > 1.5 Mpc with ^agn^dm'' where &agn is the 
AGN bias parameter and Cdm* ^^"^ halo term of the 

dark matter correlation function modeled as the Fourier 
Transform of the linear power spectrum. 

In this paper, instead, we use the HOD modeling de- 



veloped in paper II to derive the bias parameter down 
to Tp > 0.7 h^^ Mpc, instead of limiting ourselves to 
^ 1-5 Mpc. This allows for better constraints on 
6a, especially in cases where the two halo dominated (lin- 
ear) regime extends below rp « 1.5 Mpc, and a bet- 
ter treatment of the cases where the one halo term contri- 
bution is still important at rp > 1.5 /i"^ Mpc. Thus, in 
our approach, the main constraints are derived from the 
two halo term, while including the one halo term contri- 
bution in the model serves as a first-order perturbation 
from linear theory. 

Paper II discusses the HOD modeling of the CCF be- 
tween RASS/SDSS AGNs and LRGs at 0.16 <z< 0.36 
for three different models: in model A all AGNs are satel- 
lites within the same DMH as the LRGs, while models 
B and C include different realizations of the cases where 
central and satellite AGNs are included and explicitly 
parameterized. See paper II for details of these models. 
We repeat this exercise here for the AGN samples used in 
this paper to derive their bias parameters. The detailed 
results of the extensive HOD modeling will be presented 
in a separate paper (Miyaji et al. in preparation), where 
a number of new improvements in the modeling over that 
presented in paper II will be included. In this paper, we 
follow exactly the method presented in paper II. Here we 
reiterate the main procedure: 

1. Determine the central and satellite HODs of the 
tracer set galaxies ((A^g,c)(-Wii) and {NQ^s){Mh), 
respectively) from their ACF. The space density 
constraint is additionally used for volume-limited 
samples (the SDSS main galaxy and the LRG sam- 
pie). 

2. Using the derived tracer set HODs ((A^g,c)(A4) & 
(A^G,s)(Afh)) and using a parameterized model of 
the AGN central and satehite HODs ((7VA,c)(Mh) 
& {NA,s){Mh)), we fit the measured CCF between 
the AGN sample and the tracer set to constrain the 
AGN HODs. Since the galaxy ACFs have much 
higher statistical accuracy, the uncertainties in the 
tracer set HODs are negligible compared to the 
AGN HOD constraints. 

3. Derive the bias parameter of the AGN sample us- 
ing: 

^ /&h(Mh)(jVA)(Mh)0(Mh)dAfh , . 

/(A^A)(M,O0(Mh)dMh ' 

where (iVA)(Mh) = (iVA,c)(Afh) + (^A,s)(Mh), 
b\i{M\i) is the bias of DMHs with a mass Mh, and 
(/)(Mh) is the DMH mass fu nction. We use Equa- 
tion 8 of ISheth et al.l (|2001| ) with parameters from 
[Tinker et al.l (|2005[) for 6h(Mh). The la uncertain- 
ties of &A corresponds to the Ax^ < 1 region in the 
parameter space of the {Na){M-i^) model. 

5.1. The HODs of the tracer sets 

5.1.1. SDSS Mam Galaxies (0.07 <z < 0.16) 

As explained above in Section 12. 1[ the tracer set for 
the low redshift AGN samples is selected from the SDSS 
main galaxy sample with an absolute magnitude range 
of -22.1 < A/°-i < -21.1. The number density of this 
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TABLE 3 

Power, Law Fits to the Inferred AGN ACF and Derived Quantities. 



Sample 


Redshift 


Median rg 


7 


''0,7=1.9 


b{z) 


biz) 


l°g f DMH 


Name 


Range 


Zoff {h-'^ Mpc) 




(fe-i Mpc) 


PL-fits 


HOD 


Me) 



X-ray Selected AGN - RASS/SDSS AGN 



total RASS-AGN 0,07 - 0,16 

total RASS-AGN(rq) 0,07 - 0.16 

low Lx RASS-AGN 0,07 - 0.16 

high Lx RASS-AGN 0.07 - 0.16 

low Lx RASS-AGN (rq) 0.07 - 0.16 

high Lx RASS-AGN (rq) 0.07 - 0.16 

narrow line RASS-AGN 0.07-0.16 

total RASS-AGN 0.16 - 0.36 

total RASS-AGN(rq) 0.16 - 0.36 

low Lx RASS-AGN 0.16 - 0.36 

high Lx RASS-AGN 0.16 - 0.36 

low Lx RASS-AGN (rq) 0.16 - 0.36 

high Lx RASS-AGN (rq) 0.16 - 0.36 



0.13 
0.14 
0.13 
0.14 
0.12 
0.14 
0.13 

0.27 
0.26 
0.24 
0.31 
0.24 
0.30 



3.84"* 



4 96+°-*'^ 

4+0.52 
^-0.59 
c qK-l-0.65 

*-^^-0.91 
4g7+1.07 
^■"'-1,34 
o 14+O.8I 



l-^lJ-O-io 

l-»0_o.i4 

^■^'-'-0.20 
=+0.11 



1.66 
2.11 
1.61 
1 



0,12 

+0.30 
0,23 
+0.18 
0,17 
+0.41 
-0.28 



4.84 
4.33 
3.85 
4.94 
4.28 
3.86 
3.10 



+0.38 
0.41 
+0.62 
0,71 
+0.55 
0.63 
+0.56 
0.63 
+0.79 
0,95 
+0.89 
1,13 
+0.70 
0.88 



4 05+"-*^ 

1+0.40 
-0.46 

3 16+°-^l 
'-'■^"-0.70 
q+0.65 
-0.87 
1+0.63 



1 63+"-i2 
-^■"■^-0,12 



4.50"* 



5.06"* 



3 62"* 
'-'■"^-0,72 

fi OO"'"" ''^ 



1.73 
1.78 
1.82 
1.85 
1.97 



+0.12 
0.12 
+0.43 

0,28 
+0,17 
0,18 
+0.37 
0,26 
+0.18 
0.18 



4 07+0.35 
^■'^'-0.38 
-1+0.36 
-0.39 
n 10+0. 42 
'-'■-^°-0.48 
3+0.58 
-0.65 
3 5Q+0-59 

1+0.66 
-0.74 



4.50"* 



5.19"* 



5.90"* 



0,08 
0,09 
+0,15 
0,17 
+0,13 
0,15 
+0,12 
0,15 
0,19 
0,21 
+0,22 
0,26 
0,16 
0,23 

0,09 
0,11 
+0,09 
0,11 
+0,11 
0,13 
1 34+0,19 
^•■-'^-0,22 

95+°-'^'^ 
^■^^-0.17 
^+0.22 
-0,25 



1,19 
1.08 
0.98 
1.25 
1.07 
1.10 
0.80 

1.06 
1.16 
0.85 



1.63^ 



1.23 
1.25 
1.11 
1.33 
1.17 
1.36 
1.24 



+0.09 
0.08 
+0.11 
0,08 
+0.11 
0.05 
+0.16 
0.14 
+0.11 
0.08 
+0.24 
0,18 
+0.18 
0.19 



1 99+0, 13 

13 94+0-14 
^■^•^^-0.12 

1 3 ns"^"'^* 

-,0 ot:+0,17 

10.0,J_g 20 

13 I'i+O-lS 

1 00+0.25 
±o.oo_g 24 
1 3 24+0-22 



1 30+0- 10 

-^■'^'^-0-08 
7+0.08 
-0.08 
1 17+0.14 
^■^'-0.12 
S+0.12 
-0.15 
-, 10+0. 14 
l-^*-0-13 
-1+0.20 
-0-13 



1.27"* 



1.48"* 



1.50"* 



10 1 7+0-12 

I'^-l'-O-ll 

13.14 
13.01 

13 34+0-12 

13 03+0-20 
±o.uo_g 24 

1 07+O- 18 
1J-J'_0.14 



+0-11 
0-12 

+0-21 
0-22 



total RASS-AGN 





36 - 


- 


50 





42 


3.24 


+0 
-2 


96 

13 


1 


59 


+0 
-0 


40 
41 


3 


83 


+0 
-0 


55 

63 





96t« 


22 
54 


1 


02 


+0 
-0 


14 
09 


12 




28 
25 


total RASS-AGN (rq) 





36 - 


- 


50 





42 


4.13 


+0 
-1 


81 

94 


1 


98 


+0 
-0 


59 
56 


4 


01 


+0 
-0 


80 
98 


1 


20l« 


34 
47 





99 


+0 
-0 


11 

06 


12 


43l« 


25 
17 












Optically-Selected AGN 




SDSS AGN 


























total SDSS-AGN 





07- 


- 


16 





14 


3.90 


+0 
-0 


62 
70 


1 


99 


+0 
-0 


29 
25 


3 


94 


+0 
-0 


62 
72 







15 
17 





95 


+0 
-0 


17 
10 


12 




37 
31 


total SDSS-AGN(rq) 





07- 


- 


16 





15 


4.16 


+0 

-1 


85 
07 


2 


12 


+0 
-0 


55 
39 


4 


26 


+0 
-1 


83 
02 


1 


08t« 


20 
23 





94 


+0 
-0 


24 
13 


12 


63t« 


50 
45 


total SDSS-AGN 





16 - 


- 


36 





31 


4.80 


+0 
-0 


24 
27 


1 


79 


+0 
-0 


09 
09 


4 


81 


+0 
-0 


23 
24 


1 


26t« 


07 
07 


1 


29 


+0 
-0 


05 
05 


13 


lltS 


07 
07 


total SDSS-AGN(rq) 





16 - 


- 


36 





32 


4.78 


+0 

-0 


29 

33 


1 


72 


+0 

-0 


10 
10 


4 


77 


+0 
-0 


26 
27 


1 




07 
07 


1 


31 


+0 
-0 


08 
06 


13 




10 
09 


total SDSS-AGN (noX) 





16 - 


- 


36 





32 


4.95 


+0 
-0 


31 

33 


1 


85 


+0 
-0 


10 
10 


4 


94 


+0 
-0 


30 
32 


1 




09 
09 


1 


26 


+0 
-0 


08 
06 


13 


06t« 


11 
09 


total SDSS-AGN (rq+noX) 





16 - 


- 


36 





32 


4.80 


+0 

-0 


33 
35 


1 


78 


+0 

-0 


12 
11 


4 


72 


+0 
~o 


31 

33 


1 


27l« 


07 
09 


1 


29 


+0 
-0 


10 
07 


13 


lots 


13 
10 


total SDSS-AGN (onlyX) 





16 - 


- 


36 





29 


4.70 


+0 
-0 


50 
60 


1 


71 


+0 
-0 


12 
12 


4 


77 


+0 
-0 


43 
47 


1 


22t« 


12 
15 


1 


35 


+0 
-0 


10 

11 


13 




12 
15 


low Mi SDSS-AGN 





16 - 


- 


36 





31 


4.73 


+0 
-0 


36 
43 


1 


72 


+0 

-0 


14 
13 


4 


73 


+0 
-0 


32 

34 


1 


23l« 


09 
10 


1 


31 


+0 
-0 


11 

09 


13 


14lS 


13 
13 


high Mi SDSS-AGN 





16 - 


- 


36 





32 


4.36 


+0 
-0 


35 
43 


1 


75 


+0 
-0 


12 
12 


4 


48 


+0 
-0 


30 
32 


1 




08 
11 


1 


24 


+0 
-0 


08 
08 


13 


03t« 


11 
13 


low Mi SDSS-AGN (rq) 





16 - 


- 


36 





31 


4.39 


+0 
-0 


44 
60 


1 


65 


+0 
-0 


16 
15 


4 


50 


+0 

~o 


34 
37 


1 




10 
13 


1 


30 


+0 

-0 


09 
12 


13 


13lS 


11 
18 


high Mi SDSS-AGN (rq) 





16 - 


- 


36 





32 


4.53 


+0 
-0 


44 
56 


1 


70 


+0 
-0 


14 
13 


4 


64 


+0 
-0 


37 
40 


1 


19tS 


10 
13 


1 


30 


+0 
-0 


12 
08 


13 




14 
12 


total SDSS-AGN 





36 - 


- 


50 





42 


4.41 


+0 
-0 


36 

44 


1 


69 


+0 
-0 


10 
10 


4 


58 


+0 
-0 


29 

31 


1 


23l« 


09 
11 


1 


33 


+0 
-0 


07 
08 


13 




09 
12 


total SDSS-AGN (rq) 





36 - 


- 


50 





42 


4.66 


+0 
-0 


36 
43 


1 


78 


+0 
-0 


11 

11 


4 


76 


+0 
-0 


32 
34 


1 


30l« 


09 
12 


1 


32 


+0 
-0 


08 
08 


13 


03l« 


11 
11 


total SDSS-AGN (noX) 





36 - 


- 


50 





42 


2.93 


+0 
-0 


51 

66 


1 


39 


+0 
-0 


10 
10 


3 


16 


+0 
-0 


28 
30 





92l« 


09 
12 


1 


42 


+0 
-0 


08 
09 


13 




09 
11 


total SDSS-AGN (rq-|-noX) 





36 - 


- 


50 





42 


3.79 


+0 
-0 


46 

58 


1 


52 


+0 
-0 


10 
10 


3 


80 


+0 
-0 


32 
35 


1 


09t« 


09 
11 


1 


43 


+0 
-0 


10 

09 


13 




11 
11 


low Mi SDSS-AGN 





36 - 


- 


50 





40 


4.59 


+0 
-0 


47 
61 


1 


72 


+0 
-0 


16 
16 


4 


70 


+0 
-0 


40 

43 


1 


26l« 


12 
15 


1 


29 


+0 
-0 


11 

09 


13 


02l« 


14 
14 


high Mi SDSS-AGN 





36 - 


- 


50 





43 


2.92 


+0 
-1 


66 
02 


1 


48 


+0 
-0 


16 
18 


3 


42 


+0 
-0 


38 
42 





90l« 


14 
19 


1 


27 


+0 
-0 


11 
10 


12 




15 
15 


low Mi SDSS-AGN (rq) 





36 - 


- 


50 





40 


4.68 


+0 
-0 


48 
61 


1 


82 


+0 
-0 


18 
17 


4 


74 


+0 
-0 


45 
49 


1 


311^ 


14 
17 


1 


26 


+0 
-0 


13 

12 


12 




18 
19 


high Mi SDSS-AGN (rq) 





36 - 


- 


50 





44 


3.59 


+0 

~o 


61 

92 


1 


56 


+0 

-0 


16 
17 


3 


87 


+0 
-0 


42 
47 


1 


05tS 


14 
19 


1 


35 


+0 
-0 


12 

10 


13. 


05t«: 


14 
14 



Note. — Values of tq, 7, and ro,-y=i.9 are obtained from a fit to Wp{rp) over the range rp = 0.3-15 h ^ 
Mpc for all samples using the full error covariance matrix and minimizing the correlated values. The given 
bias parameters, b — (78,AG'Niz)/ag,iz), are based on the best power law fit parameter and from HOD modeling 

(HOD). To derive log M^^-^, we use the bias parameter from HOD modeling. See Table [T] for the definition of 
the samples. 



sample is (9.77 ± 0.16) x 10""^ ^V^'^^ 
five-parameter model bv lZheng et al.l |2007D to represent 
the central and satellite HODs of our low redshift tracer 
set: 



(7VG,c)(Mh) = - 



1 + erf 

(iVG.s)(Mh) = (iVG,c)(Mh) 



log Mh — log M„ 

flog M 



(8) 



This form involves a step function with a lower mass 



cutoff Mmin, which is smoothed by incorporating the er- 
ror function (erf) with the width of the cutoff profile 
Clog A/. For the detailed description o f the d ifferent pa- 
rameters, see Sect. 3.2 of iZheng et al.l ()2007f) . A limita- 
tion of our current fitting software is that it allows a max- 
imum of two simultaneous variable parameters. Thus, we 
search for acceptable fits in a two parameter space, while 
fixing other parameters to reasonable values. We note 
that, for our current purposes in this paper of obtaining 
correct 6a values, which is mainly constrained by the two 
halo term, the detail of the HOD is not critical (see dis- 
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Fig. 8. — Projected tracer set ACFs and CCFs between the different total AGN samples and tracer sets in the corresponding redshift 
ranges. The order of the panels is the same as in Fig. [7] Each ACF/CCF plot contains the corresponding best-fit HOD model (blue solid 
line) along with the one halo (red dashed line) and two halo (green long-dashed line) terms and shows the residuals (data/model) below 
each individual plot. We show the rp range that is used for the fits, i.e., rp > 0.2 Mpc and rp > 0.7 Mpc for th e tracer ACFs and 
the AGN-tracer CCFs respectively. The LRG ACF data points shown as crosses are not used in the fit (see Sect. [5~1.2| |. 

and (TiogM = 0.39. We then search for the best fit model 
for our measurement of the SDSS main galaxy ACF in 
the remaining two parameters (log Mmin and as) by min- 
imizing the correlated x^, taking into account the density 
constraint (see Eq. 18 of paper II). We find an excellent fit 
to the data with the best fit values and la uncertainties of 
log(Mmi„/[/i-^ Mq]) = 12.81±0.01 and Us = 1.16±0.02, 
with an associated bias parameter of 6g = 1-43 ± 0.01. 
The uncerta i nties i n these values are lower than those in 
iZheng et al.l (|2007f) as we perform the fit with only two 
free parameters instead of five. For the fit, we use ?x'p(rp) 



cussions below in Sect. l5.2p . A Code that searches for the 
best-fit values and confidence ranges of model parame- 
ters in a many-parameter space using the Markov-Chain 
Monte-Carlo (MCMC) is under development and will be 
used in a future paper. 

Zheng et al. (2007.) perform fits to the full five pa- 
rameter space for several luminosity-threshold subsets 
of the SDSS main galaxy sample. We take advantage 
of the results of their < 21.0 sample, which has 
a similar selection criteria to our SDSS main galaxy 
sample, and use their results to fix three parameters: 
\og{Mo/[h~^ Mq]) ^ 11.92, log(M{/[/i-iM0]) = 13.94, 



measurements in the range 0.2 < rp[h ^ Mpc] < 40. 
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First, we follow our HOD modeling described in paper II 
and exclude data points that fall in the transition region 
between the one halo dominated and two halo dominated 
regimes. In the case of the SDSS main galaxy sample, 
using scales of 0.4 < rp[h'^ Mpc] < 1 results in essen- 
tially the identical best-fit HOD. The SDSS main galaxy 
ACF and the best fit model are shown in Fig. [S] (upper 
left panel). 

5.1.2. Lununous Red Galaxies (t).16 < z < 0.36^ 

The tracer set in the intermediate redshift range is the 
LRG sample with -23.2 < M^-^ < -21.2 (Sect. [221). 
The derivation of central and satellite HODs of this sam- 
ple is described in detail in paper II . In su mmary, we 
start from the results of 'Z heng et al.l ([2009) and make 
a two parameter adjustment to find the HOD that fits 
best to our LRG ACF, including the number density 
constraint. For satellite LRGs, we interpolate between 
the Ay^ = 4 upper and lo wer bounds of the satellite 
HOD by 'Zheng^etlSI (|2009D . while for the central LRG 
HOD we shift their central HOD horizontally along the 
logAfh axis. First, we fit Wp(r-p) measurements in the 
0.2 < Tp [h^^ Mpc] < 40 range. Significant residuals re- 
main in the transition region between the one halo and 
two halo term dominated regimes, due to the fact that 
our HOD modeling neglect the effects of halo-halo colli- 
sions in the two halo term. Unlike with the main galaxy 
sample, this effect is not negligible with the LRG sample, 
due to a larger transition region. We therefore neglect 
data points in a range of 0.46 < [h~^ Mpc] < 2.8 and 
find a good fit to the data (Fig. [51 middle left panel). The 
associated LRG bias parameter is 5g — 2.20±0.01, where 
this error includes only the statistical \a uncertainty of 
the fit. 

5.1.3. Extended LRGs (1D.36 <z< 0.50; 

Unlike for the cases of the SDSS main galaxy and LRG 
samples, there is no template HOD model in the litera- 
ture for a sample with almost identical selection criteria 
to those for our extended LRG sample, i.e., a non- volume 
limited sample of LRGs with -23.2 < M°-3 < -21.7 at 
0.36 < 2 < 0.50. 

The closest sample that we can use as our template 
is the -2 3.2 < M°-^ < -2 1.8, z - 0.3 LRG sample, 
for which iZheng et al.l (|2009( ) made a detailed HOD in- 
vestigation, in addition to the -23.2 < M°-^ < -21.2 
LRG sample. Thus, we follow the approach in pa- 
per II for the LRG sample and search for the best-fit 
HOD model of the ACF o f our extended LRG sample 
by adjusting IZheng et all (|2009D HOD resuhs for the 
-23.2 < < -21.8 LRG sample. 

In short, we take the c e ntral and satellite HODs from 
Fig. 1 (b) of IZheng etHl (|2009D and search for the best- 
fit HOD by tweaking the template. We shift their cen- 
tral HOD horizontally by d in logMh and interpolate 
between their upper and lower bounds (A^^) of their 
satellite HOD, with the dividing ratio of / : (1 — /), 
where / = (/ = 1) represents their lower (upper) 
bound on the satellite HOD (see paper II for details). 
We search for the best-fit HOD model in the (/, d) space, 
where the HOD-model predicted Wp(rp) function is cal- 
culated at z = 0.42. As with the LRGs, we exclude the 
0.46 < rp[h^^ Mpc] < 2.8 range from the fit. We obtain 



the best-fit with / = 4.2 ± 0.8 and d = -0.29 ± 0.04. We 
do not use the density constraint in the fit, because the 
extended LRG sample is not volume-limited and there- 
fore the number density is not accurately determined. 
However, the best-fit model gives a density of 4.3 x 10~^ 
Mpc~'^, which is close to the number density of the 
extended LRGs calculated in the volume-limited portion 
of the sample (see Table ^ . 

5.2. AGN biases and Typical Halo Masses 

For the AGN HOD model, we use model B of paper II, 
which parameterizes the number of central ((A^a c)) and 
satellite ((A^a.s)) AGNs in a DMH: 



{Na, 



,) = fAeiAh 



Mmin), 

M„i„)(Mh/Mi)"% 



(9) 



where Q{x) is the step function (equal to 1 at a; > 0; 
at a; < 0) and /a represents the AGN fraction (duty 
cycle) among central galaxies at Mh > Mmin . We 
use log Mi/Mmin — 1.36 in Eq. [9l which Zehav i et al.l 
(|2005bf) find to be a typical value at which a DMH hosts 
on average one satellite galaxy in addition to a central 
galaxy. We model the CCFs between our tracer sets and 
the AGNs using the same method as in paper II: 
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[Pm{k) + P2i,{k)]Jo{krp)dk (12) 



where nAijiQ) is the number density of AGNs (tracers), 
y{k,M\^) is the Fourier t ransform of the Nava rro, Frenk, 
& White profile rNFW: iNarorro et al.l[l997h . Pii„(fc) is 
the linear power spe ctrum of the den s ity fie ld with the 
transfer function bv lEisenstein fc Hul ()1998() . and Jq{x) 
is the zeroth-order Bessel function of the first kind. 

Due to the low space density of AGNs, most CCFs 
do not have a signal-to-noise ratio on small scales that is 
sufficient for applying the statistics. While a majority 
of CCFs contain > 16 pairs per bin at > 0.3 Mpc, 
we have to exclude rp < 0.7 Mpc bins for several 
CCFs in order to have at least 16 AGN-galaxy pairs per 
bin. To derive 6a in a consistent way for all CCFs, we 
use only rp > 0.7 Mpc bins to derive column 'b{z) 
HOD' for aU of the AGN samples given in Table [3 

Since the purpose of using the HOD modeling in this 
paper is to derive reliable values of 6a, we do not dis- 
cuss detailed results using other models, which will pre- 
sented in Miyaji et al. (in preparation). The derived 
bias parameter, which is mainly constrained by the two 
halo term, is not very sensitive to our particular choice 
of HOD model. To verfiy this, we repeat the HOD 
fits to the CCFs with different values of the parame- 
ter log Mi /Mmin (Eq. [H). For various SPSS lum inosity- 
threshold galaxy samples, iZehavi et al.l ()2005bf) found 
the range of this parameter to be 1.0 < logMi/Mmin < 
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Fig. 9. — HOD modeling bias parameter ^agn ^ ^ function of 
redshift for our X-ray selected RASS/SDSS samples (open sym- 
bols) and optically-selected SDSS AGN samples (filled symbols). 
Black symbols represent the total AGN samples, while colored 
symbols represent subsamples (low and high Lx in green and yel- 
low, radio-quiet in blue). The dotted lines show the expected h{z) 
of typica l DMH masses Mdmh based on Sheth et al. (2001) and 
[Tinker et al., 12003) . where masses are given in log Mdmh in units 
of Mq. For visualization purposes, we slightly shift the red- 
shifts of some AGN samples. 

1.5. Thus we fix logMi/Mmin to 1.0, 1.36 (our default), 
and 1.6 and obtain 6a values in each case. The best fit 
values of &a typically vary only by « 0.01 among these 
three cases, demonstrating the robustness of deriving the 
bias parameter using this method. 

5.3. AGN HOD bias results 

In Fig. [9] we present the main results of our study. We 
show the HOD bias parameter for our different X-ray and 
optically-selected AGN samples as a function of redshift. 
All AGN samples are consistent with a host DMH mass of 
log (Mdmh/[^"^ Mq]) ~ 12.4 - 13.4 (see Tabled. Sam- 
ples in which the radio-loud AGNs have been excluded 
have very similar clustering amplitudes as the total sam- 
ples, in all three rcdshift ranges. The clustering signal 
of narrow-line RASS/SDSS AGNs at 0.07 < ^ < 0.16 is 
also very similar to broad-line RASS/SDSS AGNs at the 
same redshift. Furthermore, weak X-ray luminosity de- 
pendences on the broad-line AGN clustering amplitude 
are found at both 0.07 < z < 0.16 and 0.16 < z < 0.36. 

5.4. Power Law versus HOD Derived Bias Parameters 

The various AGN samples studied here allow us to 
compare the power law fit derived bias parameters with 
those from HOD modeling. In Fig. [10] we plot the bias 
values from the power law fits and HOD modeling listed 
in Table |3l The 6(H0D) values are often larger than the 
6(power law fit) values. Furthermore, the uncertainties 
on the HOD bias are often lower than on the power law 
bias. This is because the power law fits are based on 
the inferred ACF, while the HOD modeling is directly 
applied to the CCF which has lower statistical uncer- 
tainties. 

The relatively narrower distribution and on-average 
higher &(HOD) values compared to the &(power law fit) 
values are caused by strong variations between the sam- 
ples in the one halo term, while the variations in the two 
halo term are smaller. As described in the beginning of 
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Fig. 10. — Comparison between bias parameters derived from 
power law fit versus HOD modeling. We highlight the most im- 
portant samples at the three redshift ranges: total X-ray selected 
RASS/SDSS AGN samples (open triangles) and total optically se- 
lected SDSS AGN samples (filled triangles). The la uncertainties 
for the different methods are plotted as error bars. The dotted line 
shows a 1:1 correspondence. 



Section [5l power law fit bias measurements commonly 
use smaller scales that are in the one halo term (our fit- 
ted range is rp = 0.3 — 15 Mpc) in order to increase 
the statistical significance. If power law fits are restricted 
only to larger scales, the method suffers from the problem 
that the lowest scale where the linear biasing scheme can 
still be applied varies from sample to sample and remains 
ambiguous. HOD modeling allows, in principle, the use 
of the full range of scales since the method first deter- 
mines the one and two halo terms and then constrains 
the linear using data down to the smallest rp values that 
are dominated by the two halo term for each individual 
sample. 

Therefore, an apparently low power law bias value of 
one sample compared to another can be due to a lower 
clustering strength in the one halo term, while the "true" 
linear bias parameters may be similar. If we use scales 
of rp > 1.5 h~^Mpc for the power law fits, the errors on 
the bias are much larger, but the values are statistically 
consistent which those derived from the HOD model fits. 

In summary, on small (< Ih'^ Mpc) and large (> 
lOh-^ Mpc) scales a power law model is not a good fit to 
the data. HOD modeling is currently the optimal method 
to establish the large-scale bias parameter, provided the 
adopted HOD parameterization can adequately describe 
the true galaxy and AGN DMH population. With larger 
data sets in the future, the uncertainties in the mea- 
surements will be decreased further. To derive a reliable 
picture of AGN clustering, bias parameters should be 
inferred from HOD modeling, or at least from the com- 
parison of the correlation function with that of the dar k 
matter only in the linear regime (jAllevato et al.ll20lH) . 
because systematic errors based on power law bias pa- 
rameters will be larger than the statistical uncertainties 
of the clustering measurement. In the following discus- 
sions, we use the HOD model bias parameters whenever 
we compare our different AGN samples. In the case 
where we compare our results to other studies, which 
list only power law fit bias values, we also use power law 
fit bias values to be consistent. 
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6. DISCUSSION 

Our clustering measurements of luminous broad-line 
AGNs yield three independent data points in the poorly- 
studied low-redshift range for both X-ray and optically- 
selected broad-line AGNs. In addition, we measure the 
clustering signal of X-ray detected narrow-line AGNs. 
We derive the bias parameter of the different samples 
based on power law fits and HOD modeling. 

6.1. Comparison with X-ray selected broad-line AGN 
Clustering Measurements 

There are no published clustering measurements of X- 
ray selected broad-line AGNs with comparable low un- 
certai nties at low redshifts. For example, Mullis et al.l 
measure the clustering strength of broad-line 
AGNs in the ROSAT^'EP survey, for which we derive a 
bias parameter of h = l-SSj^^ {z = 0.22, paper I). 

At higher redshift lAllevato et al.l (|2011l ) compute the 
clustering strength of X-ray (0.5-2 keV) unabsorbed & 
absorbed and narrow & broad-line AGN in the XMM- 
Newton COSMOS field in different bins where the me- 
dian redshift of the subsamples varies from z = 0.5 
to z = 2.5. They find that broad-line AGNs reside 
in DMH of log {MY)ym/[h-^ Mq\) - 13.2, indepen- 
dent of redshift. The average luminosity of their broad- 
line AGN sample is log (L2_iokGv/[ergs~^]) ~ 44.1 
(intrinsic absorption-corrected luminosity; V. Allevato, 
private communication 2011) which corresponds to log 
(I/o.i-2.4kGv/[ergs~^]) 44.2 assuming F = 1.8. Conse- 
quent ly, the broad-line AGNs studied bv lAllevato et al.l 
(|2011l ) have X-ray luminosities comparable to our total 
RASS/SDSS AGN sample at 0.16 < z < 0.36. The de- 
rived DMH mass of log (Mdmh/[/i"^ M©]) - 13.2 for 
our low-redshift RASS/SDSS AGN sample is in excel- 
lent agreement with the v alue found at z ^ 1.5 for the 
broad-line AGN sample of lAllevato et al.l (|2011f ). 

6.2. Optical versus X-ray selected ACN clustering 
properties 

Our study allows us to compare X-ray and optically- 
selected AGN clustering measurements over three inde- 
pendent low redshift ranges. As the same procedure and 
identical tracer sets are used for calculating the cluster- 
ing strengths of the X-ray and optically-selected AGN 
samples, we expect there to be minimal systematic er- 
rors when comparing the clustering properties among our 
different AGN samples. 

Since most of the results in the literature are based on 
power law fits to the correlation function, we compare our 
power law bias parameters from various clustering mea- 
surements of X-ray and optically-selected SDSS AGN 
samples and galaxy clustering measurements in Fig. 1111 
The bias values of our X-ray and optically-selected AGN 
populations clearly fall within the region occupied by 
~L* galaxies. 

We do not find a very significant difference between the 
clustering amplitude of our broad-line optical and X-ray 
AGN samples in the different redshift ranges considered 
here. The differences in the HOD model bias parame- 
ters between the total X-ray and optically-selected AGN 
samples in our three (increasing) redshift ranges are 1.5f7, 
O.lcr, and 2. Oct. As our samples at 0.16 < z < 0.36 have 
many more X-ray and optically-selected AGNs than the 



samples at lower and higher redshifts, we consider the 
results in this redshift range to be the most reliable. We 
do not find a significant difference in the clustering of 
X-ray and optically-selected AGN in this redshift range. 
Moreover, given that only one out of three sample has 
a difference of ^2. Oct, we conclude that it is likely not 
significant. 

Although the RASS/SDSS AGN samples extend to 
lower AGN luminosities than the optically-selected SDSS 
AGNs, both samples span roughly the same luminos- 
ity range for broad-line AGNs, which may account for 
the similar clustering properties. RASS and SDSS are 
also well matched in ter ms of the depth and s e lectio n 
of broad-line AGNs (see lAnderson etldl l200l l200l . 
Optical surveys are known to find at brighter magni- 
tudes predominantly X-ray unabsorbed broad-emission 
line AGNs. ROSATa soft energy range allows primarily 
for the detection of X-ray unabsorbed AGN and, hence, 
is also biased also toward broad-line AGNs. Therefore, 
we do not expect strong systematic differences between 
these sampes due to the RASS/SDSS AGN selection. 

For the optically-selected SDSS broad-line AGNs in 
the redshift ranges of 0.16 < z < 0.36 and 0.36 < 
z < 0.50, the sample is large enough to create sub- 
samples of optically-selected AGN that are not detected 
in RASS. Furthermore, at 0.16 < z < 0.36 we select 
optically-selected AGNs that are also detected by RASS. 
We use these samples as a consistency check to verify that 
the clustering properties between X-ray and optically- 
selected AGN samples are not significantly different. As 
shown in Table 3 all of these AGN samples agree well 
within their Ict uncertainties in their HOD model bias pa- 
rameters. This provides compelling evidence that there 
is indeed no strong difference between broad-line X-ray 
and optically-selected AGNs at similar luminosities at 
low redshifts. 

Figure [TT] compares our AGN clustering results to 
other X-ray and optically-selected AGN clustering stud- 
ies. The properties of vari ous clusterin g stud ies are 
given in Table 4 of paper I. iHickox et all (|2009t ) study 
the clustering properties of AGNs in the AGES sur- 
vey. As they publish only the AGN CCFs with galax- 
ies, we use their best power law fits (R. Hickox, pri- 
vate communication 2011) to infer their AGN ACF and 
the bias parameter, following our approach described in 
Sect. 14. 1[ From th e redshift distributions presented in 
IHickox et al.l (|2009D, we compute the effective median 
redshift for the CCFs — 0.37) and derive a bias 

value of &x,Hi09 = l-20+^:ra. 

Our finding of detecting no significant difference in the 
AGN clustering properties between X-ray and optically- 
selected AGNs at low redshifts may appear to be in con- 
trast to AGN clustering measurements at higher redshifts 
(z > 0.7), where optically-selected AGN samples have a 
lower clustering strength than X-ray selected AGN sam- 
ples (Fig. [TT|) . Furthermore, the clustering of X-ray se- 
lected AGNs is roughly consistent with red galaxies at 
higher redshifts. However, some of the X-ray clustering 
studies significantly underestimate their uncertainties by 
using Poisson errors instead of jackknife errors. More- 
over, X-ray and optically-selected AGN samples at these 
redshifts select AGNs with different intrinsic properties. 
While the optical AGNs are mainly drawn from large sky 
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Fig. 11. — Power law bias parameter ^agn = o-8,AGn{z)/'^8D{z) as a function of redshift for various X-ray and optically-selected AGN 
samples as well as blue and red galaxies. We plot our power law fit bias parameters to compare different studies in a consistent manner. 
Black open symbols represent X-ray selected AGN samples, while black filled symbols represent optically-selected AGN samples. Clustering 
measurements for red and blue galaxies are shown as red and blue symbols at different redshifts. For the explanation of the dotted lines 
see Fig. [9] For visualization purposes, we slightly offset the redshifts of the X-ray and optically-selected SDSS AGN samples. 



area surveys and sample luminous predominantly broad- 
line AGNs, the X-ray selected AGN samples derive from 
very deep observations covering only a few square de- 
grees on the sky. Consequently, the X-ray selected AGNs 
have, on average, much lower luminosities. Additionally, 
the X-ray samples include absorbed AGNs, which results 
in a large fraction of narrow-line AGNs that are missed 
in the optical AGN samples at these redshifts. 

6.3. Broad-line versus Narrow-line AGNs 

The differences in the clustering signals between X- 
ray and optically-selected AGN samples at z > 0.7 can 
potentially be accounted for cither as the result of a large 
luminosity difference between the samples or the large 
fraction of X-ray absorbed, optically narrow-line AGNs 
in the X-ray AGN samples. 

To test the latter assumption, we measure the cluster- 
ing pro perties of narrow-line R ASS /SPSS AGNs classi- 
fied by [Anderson et all (|2007f) in the redshift ranges of 
0.07 < z < 0.16 and 0.16 < z < 0.36. The low num- 
ber density of narrow-line RASS/SDSS AGNs at 0.16 < 
z < 0.36 forces us to use scales of rp > 1.1 Mpc to 



apply the x^-statistics during the HOD modeling. The 
clustering strength of the narrow-line and total broad- 
line RASS/SDSS AGN samples at 0.07 <z< 0.16 agree 
well with each other (Table [3|), although the uncertainty 
for the narrow-line AGN sample is large. The power 
law fit bias parameter for the narrow-line RASS/SDSS 
AGN sample is significantly lower than the bias derived 
from the HOD model (see Fig. [9]) and would suggest a 
significantly lower clustering amplitude than the broad- 
hne RASS/SDSS AGN sample. This is caused by differ- 
ences between the samples on small scales where the one 
halo term dominates. The power law fits use these small 
scales, while the HOD modeling mainly relies on the two 
halo term to determine the large-scale clustering. 

At 0.16 < z < 0.36 the narrow-line AGN HOD bias 
parameter is too poorly constrained (b = I.OI^q'Iy) to 
allow for a detailed interpretation. While we cannot 
determine whether narrow-line RASS/SDSS AGNs clus- 
ter similarly or less than broad-line AGNs in that red- 
shift range, we can rule out the conclusion that they 
are significantly more clustered. Two important points 
are worth noting: first, as ROSAT can detect only 
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moderately X-ray absorbed AGNs (log {Nn / [cm.'"^]) < 
22), while XMM-Newton and Chandra are sensitive to 
much more absorbed (and lower luminosity) AGNs, the 
narrow-line AGNs detected by RASS may not be as com- 
mon. Second, the relative clustering strength of narrow- 
line AGNs could change with redshift if there is a dif- 
ference in how AGN activity is triggered at different 
cosmological epochs. However, other studies confirm 
that low-redshift narrow-line (radio-quiet) AGNs are not 
strongly clustered and are hosted in galaxies that do not 
diff er significantly from typical non-AG N galaxies, (e.g, 
iMandclbaum ct al. 20091 iLi et alll2006l) . 

lAllevato et al.. (2011) find that X-ray selected narrow- 
line AGNs in the XMM-Newton COSMOS field clus- 
ter slightly lower (2.3(t) than X-ray selected broad-line 
AGNs and reside in DMHs of log (AfDMH/[/i"^ Mq]) ~ 
13.0 in the redshift range z ^ 0.5 — 1.0. However, 
their narrow-line AGNs have an average intrinsic (ab- 
sorption corrected) log (L2-iokcv/[ergs^^]) — 43.1 (V. 
Allevato, private communication 2011) which is an or- 
der of magnitude lower than the average luminos ity of 
their broad-hne AGNs. When lAllevato et al.l (|2011[ ) con- 
sider only the X-ray properties of the sources to create 
X-ray absorbed and X-ray unabsorbed subsamples, in 
which both have an almost identical mean luminosity of 
log (L2_iokev/[ergs~^]) ~ 43.65, they find that X-ray 
absorbed AGNs cluster slightly less (2.6(t) than X-ray 
unabsorbed AGNs. 

To summarize, the difference between the AGN cluster- 
ing properties between X-ray and optically-selected AGN 
samples at z > 0.7 is likely not due to a strongly clus- 
tered population of narrow-line AGNs in the X-ray sam- 
ples. However, these objects do have significantly lower 
luminosities than optically-selected broad-line AGNs. 

6.4. Impact of Radio-loud Broad-line AGN on the 
Clustering Signal 

For each of the various broad-line AGN samples stud- 
ied here, we create subsamples where we have excluded 
radio-detected AGNs to study the impact of radio-loud 
AGNs on the derived clustering strength. We do not 
find any significant differences between the AGN sam- 
ples that include or exclude radio-detected AGNs. The 
HOD bias parameters for these samples agree well within 
their la uncertainties in all three redshift ranges studied 
(see Fig. However, the samples are similar as only 
approximately 10-20% of all broad-line AGNs have also 
radio (FIRST) detections, so this result may not be par- 
ticularly constraining. 

As mentioned in Section [321 the SDSS AGN target se- 
lection gives the highest priority to point sources that are 
detected in FIRST (above a certain fiux limit) without 
considering their colors. We select all AGNs {n ~ 187) 
with 0.16 < z < 0.36 which have the SDSS FIRST tar- 
get selection fiag equal to 1. Hence, these objects were 
only selected on the basis of having a significant FIRST 
radio flux and can be understood as a well-defined radio- 
selected AGN sample ((z)=0.28, (Af,)=-23.28). We com- 
pute the CCF of this AGN sample with the LRG tracer 
set. Due to the low number of AGNs in this sam- 
ple, the HOD bias parameter has large uncertainties. 
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Fig. 12. — Power law bias versus redshift (similar to Fig.|9ll, com- 
paring the clustering of differently-selected AGN samples (radio 
and non-radio samples) and galaxies. Note that the HOD model 
bias parameter for the narrow-line RASS/SDSS AGN sample is 



significant higher (b ■ 
here. 



1.24lj;;v;°) than the power law bias shown 



SDSS AGNs that have a radio-detection (including the 



187 radio-selected SDSS AGNs; (z)=0.27, (A/i)=-22.84) 
yields b — 0.97j;Q'i[g. Furthermore, we compute the CCF 
of only the RASS/SDSS AGNs with 0.16 < z< OM tha t 
are marked as radio sources in [Anderson et al.l (|2007f) . 
The sample, which contains 144 objects and has (2;)=0.25 
and log (Lo.i-2.4kev/[ergs~-^]))~ 44.36, yields a HOD 
bias parameter of 6 = 1.08^°' 3°. No constraining re- 
sults can be drawn by using only radio-detected SDSS 
AGN samples, because these samples contain too few 
objects. However, all values are consistent with the clus- 
tering strengths of the other AGN samples in the same 
redshift range. 

Many previous studies fe.g.. iMagliocchetti et al.ll2004l : 
iHickox et cd]|2009t [Mandelbaum et al.ll2009[ ) found that 
radio-loud AGNs cluster more strongly than AGNs with- 
out the presence of radio emission. At first glance 
our results may appear to be in contrast to these find- 
ings. However, those studies focus on the clustering 
properties of optical narrow-line (instead of broad-line) 
AGNs based on diagnostic emission-line ratios and ra- 
dio luminosity. Moreover, when comparing the cluster- 
ing strength of different samples it is essential to take 
into account the involved luminosities. Figure [121 shows 
the bias of various AGN and galaxy samples at lower 
redshift (z < 0.6), focusing on comparing radio-selected 
AGN s amples. For t h e radi o (non-broad-line) AGN sam- 
ple of IHickox et all (poM ). we derive a bias value of 
^radio,Hi09 = 2.07tQ;i3 = O'^^) by following the de- 
scription given in Sect . 6.21 

MagUocchetti et al.l ([2004) find that radio (FIRST) 
galaxies with AGN activity cluster more strongly than 
both radio-detected 2dFGRS galaxies without AGN ac- 
tivity and 2dFGRS galaxies without radio emission. It 
is not clear, however, whether the different galaxy sam- 
ples have similar luminosities and colors. If the galaxy 
samples that are used for comparison are mainly blue, 
star-forming galaxies, one would expect a significant dif- 
ference in the clustering strength based on the different 
host galaxy properties. Consequently, the clustering dif- 
ference could reflect not by AGN activity but host galaxy 
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type. The high bias parameter for radio-loud AGNs 
(6 = 2.07t[J;i3) and m oderate bias (& = 1 .20;^;^j^) for 
X-ray selected AGNs in IHickox et al.l (|2009( ) can also be 
explained by the different host galaxy populations (see 
t heir Fig. 9). 

iMandelbaum et al.l (|2009[ ) compare the clustering 
strength of optically-selected narrow-line AGNs (0.01 < 
z < 0.3 and -23 < M°-^ < -17) with and without ra- 
dio emission at /iierf stellar mass, i.e., at fixed luminosity 
and color. They find that at fixed stellar mass radio-loud 
narrow-line AGNs cluster more strongly and hence reside 
in more massive DMHs than narrow-line AGNs without 
radio emission and galaxies without the presence of AGN 
acti vity fsee Fig. [T2l). N either Magliocchctti et al. (2004) 
nor IMandelbaum et al.l (^2009) detect a difference in the 
clustering properties within the radio-loud narrow-line 
A GN sample as a fu nction of radio luminosity. 

IShen et all (|2009h co mpute the ACF of op tical SDSS 
DR5 AGNs (from the iSchneider et al.l [20071 sample) in 
the redshift range 0.4 ^ z < 2.5. As we do here, they di- 
vide their samples into FIRST-detected and undetected 
sources. Their fits (over the range 3-115 Mpc) suf- 
fer large uncertainties; therefore, they fix the slope of 
the power law fit for all samples to 7 = 2. They compute 
the ACF of the FIRST-undetected SDSS AGNs. When 
they compare this result to the CCF between FIRST- 
detected and undetected SDSS AGNs, they find a 2.5ct 
difference in vq , using only the diagonal elements of the 
covariance matrix. Fixing the slope and not using the full 
covariance matrix likely increases the systematic uncer- 
tainties in their results. By contrast, our measurements 
test a much narrower redshift range, use the full covari- 
ance matrix, allow free tq and 7 values for the po wer law 
fits, a nd use the HOD modeling approach. She n et al.l 
(|2009f ) test much brighter AGN samples at higher red- 
shifts than our samples here. The fact that they find a 
difference in the clustering of radio-detected versus non- 
detected AGNs and we do not may be a consequence of 
the samples having different luminosities or redshifts. 

If broad-line AGNs at z < 0.5 indeed have no clus- 
tering dependence on radio emission, this may suggest 
that different physical mechanisms trigger radio emis- 
sion in broad-line and narrow-line AGNs, as narrow-line 
AG Ns do have a clu stering dependence with radio emis- 
sion (|Mandelbaum et al. 2009) at fixed stellar mass. This 
would be somewhat surprising as significant radio emis- 
sion in al l AGNs is believed to be r elated to jet phenom- 
ena (e.g.. lBlandford fc Pavnd[l982l ). 

6.5. Luminosity Dependence of the Clustering Signal 

In papers I and II, we reported a possible X-ray 
luminosity-dependence in the AGN clustering strength of 
RASS/SDSS AGN in the redshift range 0.16 <z < 0.36 
(see also Fig. [SI green and yellow green symbols). The ex- 
clusion of radio-detected RASS/SDSS AGNs performed 
here tests whether the weak Lx-dependence of AGN clus- 
tering observed is due to the presence of radio-loud AGNs 
in the higher Lx sample (see Sect. I3.1TT1) . 

As shown in paper II, this luminosity-dependent clus- 
tering is more significant in the one halo term than in 
the two halo term. Using only larger scales of rp > 0.7 
Mpc for the HOD modeling in this paper (in order 
to derive bias parameter in a consistent way for all AGN 



samples) increases the error on the individual measure- 
ments. This results in a decrease in the significance of the 
luminosity-dependent clustering. We find that the 1.5a 
(2.0(7 using the power law bias of this paper) clustering 
difference between the lower and higher Lx AGN samples 
remains roughly constant when we exclude radio-detect 
AGNs from both subsamples (1.7cr using the HOD bias; 
2.4cr using the power law bias of this paper). 

The difference in the clustering strength of the lower 
and higher Lx RASS/SDSS AGN samples in the redshift 
range of 0.07 < z < 0.16 is 1.2cr (HOD model). The 
error bars on the bias of these samples are high, due 
to the low number of objects. The significance of the 
individual redshift measurements do not provide strong 
evidence for an X-ray luminosity dependence of the AGN 
clustering strength. 

Other studies fin d similar weak trends at low redshifts 
(Fig. [12]) . Ca ppelluti et al.l (12010.) compute the cluster- 
ing signal for Swift/BAT 15-55 keV selected AGNs and 
find a clustering difference of 1.6(t for the higher Lx 
AGNs relative to the lower Lx AGNs. Comparing the 
photometric galaxy density around spectroscopic SDSS 
AGNs, Serber et al. (20QE) and [Strand et al. (2008) find 
that higher luminosity AGNs have more overdense galaxy 
environments compared to lower luminosity AGNs at 
scales smaller than 0.5 Mpc. These studies are not di- 
rectly measuring AGN clustering, rather they focus on 
the immediate environments of AGNs. 

We do not detect an optical luminosity-dependence of 
the RASS/SDSS AGN clustering. Given the large sam- 
ple size and low resulting errors in these samples, this 
is a constraining result. Previous clustering measure- 
ments of optically-select ed AGNs u s ing th e 2dF QSO 
Redshift Survey (2QZ, 'Bovle et al.' '2000^, 2dF-SDSS 
LRG and QSO (2SLAQ, Cannon et al. 200_fi), and SDSS 
also find little evidence for an optical lumi nosity depen- 
dence of the AGN clustering strength fe.g.. lCroom et al] 

2002t Ida Angela et al.l 120081: iMountrichas et al.l [20091) . 
Mountrichas et al.l (12009^ use CCFs between AGNs and 
LRGs and find some indication that bright SDSS AGNs 
cluster less than faint 2SLAQ QSOs, although the re- 
sult is only margi n ally si gnificant (l.Qa). On the other 
hand, IShen et all (|2009f ) detect a stronger clustering 
strength for the 10% most luminous SDSS DR5 AGNs 
(0.4 ^ z < 2.5) at the ~2cr level. However, they caution 
that the dynamical range in luminosity probed is narrow 
and the sample size in the luminosity subsamples not 
large e nough to yield constrainin g clustering measure- 
ments. iPorciani fc NorbergI (|2006[ ) suggest that a lumi- 
nosity dependence of the clustering may be more evident 
at z > 1.3, while Ida Angela et all ()2008[ ) find hints that 
the lower redshift ranges {z < 1.3) may show more de- 
pendence with luminosity. 

The optical and X-ray luminosities of broad-line AGNs 
are connected via the optical-to-X-ray spectral index, 
which measures the ratio of the rest-frame luminosity 
density at 2500 A to 2 keV. Although the relation has 
some scatter (see [Anderson et al.l 120071 for the optical- 
to-X-ray spectral index for RASS/SDSS AGNs), more 
X-ray luminous AGNs are, on average, also intrinsically 
brighter in the optical. If indeed there is a weak X-ray lu- 
minosity dependence of the AGN clustering strength, the 
very narrow Ali range that the optical SDSS AGNs span 
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would hamper a detection. Identifying AGN activity us- 
ing X-ray emission allows us to identify low-luminosity 
AGNs, where the optical light is dominated by the host 
galaxy. This results in a wider luminosity range in X- 
ray emission than in the optical (Fig. |6l lower panel) in 
the two lowest redshift ranges. At 0.16 < z < 0.36, the 
mean optical luminosity of the high Mj SDSS AGN sam- 
ple is only a factor 2 higher (0.75 mag) than the low 
Mi SDSS AGN sample, while the X-ray luminosities be- 
tween the high and low Lx samples differ by a factor of 
4.4. This is also seen when considering the luminosity 
difference between the 90th percentile in the high lumi- 
nosity samples and the 10th percentile in the low lumi- 
nosity samples (factor / — 4.3 for optical SDSS AGNs, 
/ = 17.8 for RASS/SDSS AGNs) and covered luminos- 
ity range (/ = 30 for optical SDSS AGNs, / = 310 for 
RASS/SDSS AGNs). 

Hence, to detect a possible optical luminosity depen- 
dence in the broad-line AGN clustering strength, as 
might be expected if there is a weak X-ray luminosity 
dependence, considering the optical-to-X-ray luminos- 
ity relation, a wider optical luminosity range has to be 
tested. As these samples already include the brightest 
objects, this is only possible if one can include lower 
optical luminosities where broad-line AGNs are not ef- 
fectively selected due to an increase in the host galaxy 
starlight fraction. Only much deeper surveys at larger 
redshifts may yield the dynamical range to test the 
luminosity dependence for optically selected broad-line 
AGNs. 

The possible X-ray luminosity dependence of broad- 
line AGN clustering detected at low redshifts (in that 
more X-ray luminous AGNs are more clustered than 
less X-ray luminous AGNs) may be difficult to reconcile 
with the result that, on average, low luminosity (mainly 
narrow-line) X-ray selected AGNs at higher redshifts are 
more clustered than luminous optically-selected broad- 
line AGNs (see Fig. [Tl]) . If the X-ray luminosity depen- 
dence is real, this may suggest that different physical 
processes trigger AGN activity at different cosmological 
times or at different luminosities. 

6.6. Current Picture of AGN Clustering 

Although AGN clustering measurements are currently 
not as constraining and the interpretation of the results is 
not as clear as for galaxy clustering measurements, some 
general findings have emerged in the last few years. 

At low redshifts (z < 0.5), broad and narrow-line 
AGN cluster similarly to inactive galaxies, occupying 
DMH masses of log (Mdmh/[/i"^ M©]) ^ 12-0 - 13.5. 
This DMH mass range includes cases where the DMH is 
dominated by one >L* galaxy or a small galaxy group 
comp osed of multiple such galaxies (|Zehavi et al.ll2005bL 

Independent of the selection method, the cluster- 
ing strength of broad-line AGNs does not significantly 
change, while narrow-line AGNs show a significant in- 
crease in the clustering amplitude when radio-selected 
narrow-line AGNs are studied. Finally, more X-ray lu- 
minous broad-line AGNs may cluster more strongly than 
their lower luminosity counterparts. Although the vari- 
ous AGN samples have different luminosities, radio-loud, 
optically-selected narrow-line AGNs, very luminous X- 
ray AGNs, and red galaxies reside in somewhat similar 



high DMH masses. Lower luminosity X-ray AGNs, opti- 
cal narrow-line AGNs with no radio emission, and blue 
galaxies tend to be found in lower DMH masses. 

At high redshifts (z > 0.7), X-ray selected AGN 
samples appear to cluster more strongly than optically- 
selected AGNs. The reason for this remains unclear. 
Possibly either low-luminosity or narrow-line (X-ray ab- 
sorbed) AGNs cluster more strongly than very lumi- 
nous broad-line optical AGNs. Additionally, as some of 
the X-ray clustering studies significantly underestimate 
their systematic uncertainties it may turn out that these 
measurements are consistent with optical AGN cluster- 
ing measurements. More high-z AGN clustering mea- 
surements based on larger samples are needed to gain a 
clearer picture. 

In this paper we use AGN samples based on SDSS. 
In the very near future there is not another planned 
survey that includes photometry and a dedicated ex- 
tensive spectroscopic follow-up program for AGNs over 
such large areas as that covered by SDSS. As both large 
co-moving volumes and spectroscopic redshifts are es- 
sential for precise AGN clustering measurements, ma- 
jor improvement in our X-ray and optically-selected low 
redshift AGNs are therefore not expecte d in the very 
near future. BOSS (lEisCTistein et al.l l20Tl and BigBOSS 
(jSchlegel et all 1201 1[ ) will detect high redshift AGNs at 
z ~ 2.2, which will improve AGN clustering measure- 
ments at higher redshifts. 

In the coming years, the i? 05*^47 successor eROSITA 
(IPredehl et all 120071) wih perform an all-sky survey in 
the hard and soft X-rays, probing much fainter than 
RASS, which is expected to detect up to ~3 million 
AGNs. Additi onally, the Large Synoptic Survey Tele- 
scope (LSST, Ivezic et al.l 120081 ) is expected to identify 
-2 million AGNs in optical bands. eROSITA and LSST 
have the potential to significantly improve AGN cluster- 
ing measurements at low and high redshifts, though only 
if there are dedicated large spectroscopic follow-up pro- 
grams. 

7. CONCLUSIONS 

This work presents AGN clustering measurements at 
low redshifts in three independent redshift ranges: 0.07 < 
z < 0.16, 0.16 < z < 0.36, and 0.36 < z < 0.50. 
Extending the use of the cross-correlation method of 
iKrumpe et all (|2010a( ). we infer the auto-correlation for 
both ROSAT All-Sky Survey (RASS) and optically- 
selected SDSS broad-line AGNs. As tracer sets we use 
SDSS main galaxies, SDSS luminous red galaxies, and 
very luminous red galaxies (exte nded LRG sample). We 
apply the HOD model method (|Mivaii et all 1201 ID di- 
rectly to the measured CCFs to derive the bias parame- 
ter. We study the impact of different AGN selections on 
the clustering signal of broad- line AGNs, i.e., by exclud- 
ing radio-detected AGNs. Furthermore, we compute the 
clustering strength for RASS-selected narrow-line AGNs. 

We find no statistically convincing difference in the 
clustering of X-ray and optically-selected broad-line 
SDSS AGNs at low redshifts (z < 0.5). Diflferent AGN 
selections based on either X-ray (RASS), optical (SDSS), 
or radio (FIRST), and combinations of these, do not 
significantly change the clustering signal for broad-line 
AGNs. This appears to be in contrast to other studies 
that find stronger clustering for radio-loud AGNs (e.g.. 
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iMandelbaum et all l2009f ) . However, those results are 
based on narrow-line, low luminosity AGNs, while our 
sample consists of more luminous broad-line AGNs. For 
the X-ray selected broad-line RASS/SDSS AGNs we find 
HOD bias values of 1.23t^-^^, 1.30to^^, and 1.02i|] J* in 
the redshift ranges 0.07 < z < 0.16, 0.16 < z < 0.36, 
and 0.36 < z < 0.50, respectively, while the HOD bias 
values for the optically selected broad-line SDSS AGNs 

are 0.95to.\l, l-SSto.os, and l.3Sta°ol, respectively. The 
corresponding inferred typical dark matter halo masses 
hosting our broad-line AGNs are in the range of log 
(AfDMH/[/j~^ Mq]) ^ 12.5 — 13.2 and are consistent with 
those occupied by >L* galaxies at these redshifts. 

We measure the clustering of RASS selected narrow- 
line AGNs, which consists of a mix of NLSls, Seyfert 
1.5, 1.8, 1.9, and Seyfert 2 candidates. We do not find a 
significantly lower clustering amplitude of RASS narrow- 
line AGNs compared to broad-line AGNs, although these 
measurements are subject to large uncertainties. In ad- 
dition, we rule out that RASS narrow-line AGNs cluster 
significantly more strongly than broad-line AGNs at low 
redshifts. 

We show that the exclusion of radio-detected 
RASS/SDSS AGNs in 0.16 < z < 0.36 does not change 
the weak X-ray luminosity dependence of the AGN clus- 
tering strength that we find in paper I (in that higher 
Lx AGNs cluster more strongly than lower Lx AGNs 
at ^ 2a). We do not detect an optical luminosity de- 
pendence of the broad-line AGN clustering in the same 
redshift range, though this result is not particularly con- 
straining due to the narrow Mi range that is covered. 

We derive the bias parameter based on the best power 
law fit, the standard method used in literature, as well as 
by using HOD modeling. Important differences between 
the two techniques are found for some AGN samples. In 
particular, using a power law fit can underestimate the 
bias compared to HOD modeling. We show that HOD 
model bias parameters are more reliable and more accu- 
rately refiect the large-scale clustering strength. Larger 
AGN samples will be provided by future missions such 



as eROSITA. As these samples will have lower statisti- 
cal uncertainties, HOD model bias parameters should be 
used to avoid introducting systematic errors that could 
exceed the statistical errors and thus possibly lead to a 
misinterpretation of clustering measurements. 
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